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Abstract  of  Dissertation  Presented  to  the  Graduate  School 
of  the  University  of  Florida  in  Partial  Fulfillment  of  the 
Requirements  for  the  Degree  of  Doctor  of  Philosophy 


PROPER  MOTIONS  OF  STARS  IN  THE 
REGION  OF  THE  ORION  ASSOCIATION 

By 

Richard  Smart 
December  1993 


Chairman:  Dr.  H.  K.  Eichhorn 
Major  Department:  Astronomy 

We  present  results  from  a proper  motion  study  of  approximately  2000  stars  in  the 
Region  of  the  Orion  Association  la  and  Ib  subgroups.  Past  studies  in  this  region  have 
been  restricted  for  a number  of  reasons:  the  large  angular  size  of  the  association,  the  lack 
of  computational  power  and  the  high  observational  priority  that  the  forming  subgroups 
Ic  and  Id  have  received. 

Our  observational  material  spans  100  years,  consisting  of  photographic  plates  taken 
in  1955-6  and  1992  at  the  McCormick  Observatory  and  newly  reduced  Astrographic 
Catalogue  positions  from  around  1900.  The  second  and  third  epoch  plates  were 
measured  on  the  University  of  Minnesota  Automated  Plate  Scanner  in  two  directions 
with  two  threshold  levels.  This  provided  us  with  an  average  of  three  position  estimates 
per  image  with  an  image  centering  precision  of  1-2  microns. 

We  reduced  the  measured  rectangular  coordinates  to  equatorial  coordinates  by  the 
overlapping  plate  technique.  This  technique  takes  advantage  of  the  extensive  multiple 
coverage,  and  uses  additional  information  on  magnitude  effects  obtained  from  objective 
grating  images. 


This  study  will  provide  observations  that  can  be  used  to  restrict  the  theories  of 
stellar  association  formation,  runaway  stars,  evolutionary  ages,  pre-stellar  gas  dynamics 
and  galactic  kinematics.  Through  this  work  we  illustrate  the  potential  power  of  the 
overlapping  plate  technique  and  its  application  to  wide  field  astrometry. 
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CHAPTER  1 
INTRODUCTION 


This  study  will  make  two  major  contributions  to  astronomy:  act  as  a working 
example  of  the  power  of  using  an  overlapping  plate  reduction  technique  and  find  the 
proper  motions  of  stars  in  the  region  of  the  Orion  association.  These  observations  can 
be  used  to  increase  our  understanding  of  the  dynamics  of  stellar  associations.  This 
chapter  will  introduce  the  problem,  examine  the  previous  work,  and  briefly  outline  the 
proposed  reduction  technique  and  the  rationale  for  using  it. 

Project  Overview 

I-Orion,  at  a distance  of  450  pc  and  with  over  50  O and  B type  stars,  is  a particularly 
nearby  and  rich  stellar  association.  It  is  centered  on  B1950  right  ascension  5*’30^  and 
declination  -1°  with  a total  area  of  approximately  100  square  degrees.  The  largest 
projected  linear  diameter  of  the  association  is  about  100  pc,  and  it  has  an  estimated 
total  mass  of  about  7.6x10^  solar  masses. 

Many  authors  (cf.  Blaauw,  1964)  have  noted  that  the  study  of  association  kinematics 
provides  insight  into  many  problems  in  astronomy:  stellar  evolution,  galactic  kinematics, 
pre-stellar  gas  dynamics,  origins  of  runaway  stars  and  the  search  for  planetary  disks 
(especially  in  young  associations  / clusters).  Even  though  the  potential  for  application 
in  astronomy  is  rich,  few  concerted  studies  of  the  proper  motions  in  I-Orion  have  been 
carried  out. 

The  reason  for  this  lack  is  probably  twofold:  First,  I-Orion  is  split  into  four 
subgroups,  a-d,  and  within  these  subgroups  the  smallest  one,  Id-Orion,  is  prioritized 
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for  telescope  and  research  time  because  this  subgroup  contains  very  young  stars  that 
have  just  left  their  formative  stages  and  are  thus  important  for  the  study  of  stellar 
formation  and  the  search  for  planetary  disks.  Second,  and  probably  more  to  the  point, 
the  large  angular  size  of  the  I-Orion  association  makes  it  extremely  difficult  even  to 
obtain  accurate  relative  positions  across  the  entire  field. 

Previous  investigations  of  the  large  subgroups,  a and  b,  in  the  association  were 
therefore,  with  two  exceptions,  restricted  to  spectroscopic  and  radial  velocity  studies. 
The  exceptions  looked  at  proper  motions  found  from  different  catalogue  sources  that 
defined  the  stellar  positions  at  different  epochs.  The  inherent  nonhomogeneity  of  this 
data  — due  to  its  use  of  many  different  sources  with  systematic  and  accidental  errors  — 
has  limited  the  success  in  determining  proper  motion,  membership,  and  kinematic  ages. 
Indeed,  as  discussed  later,  the  two  most  recent  astronomical  studies  on  these  subgroups 
have  called  for  a more  accurate  proper  motion  study. 

The  present  investigation  will  correct  the  current  lack  of  accurate  proper  motion 
studies.  We  will  determine  the  proper  motions  of  stars  in  the  region  outlined  on  figure 
1 almost  complete  down  to  12th  magnitude.  To  attain  consistency  across  the  region,  we 
will  use  a sophisticated  overlapping  plate  reduction  technique  that  reduces  the  whole 
region  in  order  to  provide  positions  for  all  the  stars  simultaneously,  consistent  in  a least 
squares  sense,  at  one  epoch. 

This  is  only  possible  now  because  of  the  recent  rapid  growth  in  fast,  large  memory 
computers.  The  region  under  study  covers  part  of  both  subgroup  a and  b,  extends  over 
30  square  degrees  and  includes  more  than  2000  stars. 
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Previous  Research  on  Associations;  Rationale  for  this  Investigation 


Ever  since  the  first  determination  of  spectral  classifications  for  bright  stars,  the 
existence  of  loose  groups  of  O and  B types  stars  has  been  known  (Eddington,  1914).  In 
order  of  increasing  distance  the  groups  IC  2602,  Scorpio-Centaurus,  Il-Perseus,  I-Orion, 
I— Lacerta,  NGC  7160,  NGC  2264, 111-Cepheus,  I-Cepheus,  NGC  1502,  NGC  2169  are 
all  within  one  kiloparsec  of  the  solar  neighborhood.  Most  of  these  groups  stand  out 
conspicuously  as  bright  configurations  in  the  night  sky. 

Astronomers  carried  out  considerable  research  to  explain  the  existence  of  these 
groups.  Kapteyn  (1914)  studied  them  when  stars  were  still  believed  to  bum  helium. 
He  produced  a whole-sky  map  of  their  positions.  Collinder  (1931)  originally  classified 
subgroups  of  these  groups  as  clusters.  Dynamical  studies  by  Bok  (1934)  and  Mineur 
(1939)  showed  these  groups  must  be  unstable.  Therefore,  by  implication  the  groups 
must  be  young. 

Ambartsumian  is  credited  with  first  coining  the  name  ‘O  associations’  for  these 
groups.  In  1952  he  described  them  as  follows: 


O-associations  are  stellar  systems,  where  the  partial  density  of  B2  stars  is  larger  than  the 
average  field  density  of  these  stars,  in  such  a way  that  this  difference  cannot  be  explained  by 
chance  fluctuations;  moreover,  O or  BO  stars  are  present.  The  properties  of  O-associations  may 
be  described  as  follows: 

1.  the  linear  diameters  range  between  30  and  200  parsecs; 

2.  the  associations  contain  an  open  star  cluster  of  type  O as  nucleus; 

3.  they  include,  besides  O — B2  stars,  also  stars  of  types  later  than  B2,  sometimes  even 
Wolf-Rayet  stars,  though  it  is  difficult  to  ascertain  the  number  of  faint  stars; 

4.  sometimes  multiple  systems  of  Trapezium  type  and  star  chains  may  be  part  of  the  nuclei; 

5.  hot  giants  occur  also  outside  the  nuclei; 

6.  there  are  reasons  for  presuming  the  O-associations  to  be  unstable  systems. 
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Ambartsumian  made  extensive  investigations  of  the  properties  of  these  associations 
(cf.  Ambartsumian,  1955).  He  pointed  out  that  the  combination  of  the  work  of 
Bok,  which  showed  that  associations  are  not  bound  by  their  gravitation,  and  that 
associations  are  roughly  spherical  in  shape,  implies  that  associations  must  be  in  a state 
of  expansion.  He  estimated  a typical  association  age  to  be  ~ 10^  years.  This  kinematic 
age  agrees  well  with  the  evolutionary  age  of  the  constituent  O and  B stars.  This  work 
established  the  existence  of  O associations  as  expanding  young  groups  of  stars.  In 
1952  Markarian  catalogued  O associations  and  created  the  International  Astronomical 
Union-recommended  nomenclature  we  use  today. 

The  first  direct  observational  evidence  supporting  Ambartsumian’s  hypothesis  was 
by  Blaauw  (1952),  and  it  was  later  confirmed  by  Delhaye  and  Blaauw  (1953).  In  these 
studies,  the  authors  used  proper  motions  from  meridian-circle  catalogues  and  showed 
that  the  ( Persei  association  (Il-Perseus)  was  expanding.  The  expansion  age  found  by 
this  study  was  1.3  million  years,  again  agreeing  well  with  the  members’  evolutionary 
ages.  Blaauw  (1964)  pointed  out  that  the  systematic  errors  in  the  catalogues  used  could 
mimic  the  effect  of  expansion,  or  even  contraction.  Therefore,  the  use  of  catalogues  for 
measuring  an  association’s  expansion  is  severely  limited,  because  of  systematic  errors 
and  by  the  general  low  sky  density. 

Fredrick  (1956)  carried  out  another  investigation  of  Il-Perseus  using  McCormick 
and  Sproul  photographic  plates.  He  reduced  the  plates  using  a linear  single  plate 
reduction  with  an  averaging  of  the  measured  rectangular  positions  from  all  the  plates 
exposed  at  one  epoch.  The  total  time  base  for  the  study  was  38  years.  Although  the 
reduction  was  the  best  that  could  be  performed  with  the  theoretical  and  computational 
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tcx)ls  available  at  the  time,  his  choice  of  a simple  model  and  method  is  prone  to  many 
systematic  errors.  In  using  a linear  plate  model,  Fredrick  ignored  radial  terms  such  as 
magnitude  and  coma.  These  terms,  if  they  physically  exist  but  are  not  represented  in 
the  modelling,  can  produce  an  expansion  of  the  association  even  if  one  does  not  exist. 
However,  the  result  was  in  reasonable  agreement  with  Blauuw’s  meridian  catalogue 
study  and  Il-Perseus  was  confirmed  to  be  in  a state  of  expansion. 

This  cemented  the  observational  evidence  for  Ambartsumian’s  hypothesis  of  young, 
expanding  associations.  The  ‘discovery’  of  associations  and  their  expansion  stimulated 
research  on  these  groups.  Because  of  their  relative  youth  and  apparent  dynamical 
parameters,  studies  of  them  have  applications  to  many  other  areas  of  astronomy. 

Ambartsumian  proposed  that  associations  may  contain  subgroups.  Blaauw’s  1953 
study  on  Il-Perseus  supported  this.  In  1964  Hardie  et  al.  examined  the  spectroscopic 
properties  of  I-Orion,  the  results  of  this  study  also  suggested  subgroups.  Blauuw  in 
his  1964  review  paper  outlined  the  various  subgroups  of  local  associations  and  listed 
some  basic  properties  of  them. 

The  I-Orion  association  appeared  to  have  4 subgroups;  figure  1 shows  the  a and  b 
subgroup  boundaries.  With  a few  exceptions,  most  studies  of  Orion  have  been  limited 
to  the  subgroups  c and  d (for  a review  of  this  research  cf.  van  Altena  et  al.  , 1988).  We 
will  be  looking  at  the  subgroups  a and  b,  which  have  received  very  few  proper  motions 
studies  (a  review  of  the  literature  shows  only  two  studies:  Lesh,  1968,  and  Gieseking, 
1983).  The  large  angular  size  the  subgroups  subtend  would  have  made  an  accurate 
astrometric  reduction  at  best  difficult,  and  perhaps  impossible,  before  the  advent  of  fast, 
large  memory  computers. 
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There  have  been  some  studies  that  utilized  existing  catalogues  much  in  the  same  way 
Blaauw  carried  out  his  original  Il-Perseus  study.  One,  that  specifically  addressed  the 
la-Orion  subgroup,  was  by  Lesh  (1968).  This  study,  like  Blaauw’s,  used  meridian-circle 
catalogue  data  and  found  an  expansion  age  of  4.5  x 10^  +2.3  x 10^  years.  Primarily 
the  large  size  of  this  error  results  from  the  large  errors  in  the  proper  motions.  In 
addition,  since  the  total  number  of  stars  used  to  derive  this  age  was  only  16,  the  large 
error  also  reflects  a problem  with  small  sample  size.  In  this  study  of  I-Orion,  we  will 
determine  uniformly  derived  proper  motions  of  over  2000  stars.  The  limiting  factor  on 
the  precision  of  the  expansion  age  will  be  the  contamination  by  non-members. 

Another  study  that  incorporated  the  proper  motions  of  both  the  a and  b subgroups 
was  a membership  study  of  the  whole  Orion  association  by  Warren  and  Hessler  (1977, 
1978).  Their  study  combined  all  previous  photometric,  radial  velocity,  and  proper 
motion  measurements  and  added  new  photometric  data  to  construct  one  homogenous 
set  of  data.  They  found  expansion  ages  in  good  agreement  with  Lesh’s  and  posed 
many  unresolved  questions:  What  are  the  defining  boundaries  of  the  subgroups?  What 
are  the  distances  of  these  subgroups?  How  homogenous  is  the  association?  One  of 
the  conclusions  in  this  study  was  that  more  complete  studies  of  proper  motions 
and  radial  velocities  are  needed  for  the  association,  especially  In  the  regions  of 
subgroup  a and  b. 

Gieseking  (1983)  carried  out  an  extensive  radial  velocity  study  of  la-Orion.  He 
concluded  that  it  consists  of  two  groups,  one  of  which  may  be  rotating.  Further  analysis 
of  Hardie’s  spectroscopic  data  supported  the  two-group  hypothesis.  Like  Warren  and 
Hessler,  Giesking  also  concluded  that  “it  would  be  interesting  to  investigate  whether 
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these  kinematically  distinguished  (via  radial  velocities)  groups  may  be  recognized  also 
by  their  proper  motions  and  available  space  motions.  . . . Unfortunately  the  proper 
motions  available  are  by  no  means  adequate.”  This  study  will  provide  the  proper  motions 
for  that  investigation. 

Advances  in  Determining  Proper  Motions  and  Positions 


Proper  motions  are  the  time  derivatives  of  the  position  components  of  a star  in  an 
inertial  reference  frame  or  a frame  that  is  a good  approximation  to  one.  The  standard 
method  for  determining  proper  motions  is  to  find  the  position  of  a particular  star  at  two 
widely  spaced  epochs  in  the  same  frame  of  reference.  The  proper  motion  is  closely 
approximated  by  the  difference  of  these  positions  divided  by  the  epoch  difference.  As 
we  have  seen,  previous  methods  used  catalogues  to  define  the  stellar  positions  at  two 
epochs.  However,  this  method  will  only  provide  proper  motions  of  limited  precision. 

A more  precise  method  is  to  photograph  the  sky  at  at  least  two  epochs  and  compare 
the  positions  of  the  stars.  This  is  essentially  the  method  used  by  Fredrick  in  his  II- 
Perseus  study.  However,  as  we  have  pointed  out,  the  model  his  study  used  to  find 
equatorial  coordinates  in  the  sky’s  frame  of  reference  from  the  image’s  measured 
rectangular  coordinate  was  too  simple  to  account  for  the  possible  effects  of  a star’s 
magnitude  or  the  telescope’s  coma  on  the  star’s  rectangular  position.  (Work  carried 
out  by  Eichhom,  1956,  revealed  the  presence  of  a large  coma  term  in  the  McCormick 
telescope.)  He  made  this  simplification  because  the  numerical  manipulation  required 
to  model  magnitude  and  coma  effects  would  have  resulted  in  too  large  a problem  to 
handle  at  that  time. 
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In  addition,  the  Orion  region  in  this  study  covers  an  area  of  approximately  10 
McCormick  photographic  plates.  The  reduction  of  a number  of  non-concentric  plates 
covering  a large  area  will  lead  to  varying  degrees  of  precision  across  the  region. 
Eichhorn  (1960,  1984)  developed  an  improved  method  called  the  overlapping  plate 
technique.  In  this  method  the  observer  exposes  the  region  such  that  all  plates  overlap 
with  at  least  one  other  plate  (and  usually  more).  This  method  assures  that  there  will  be 
stars  in  common  to  those  plates;  and  because  these  plates  were  exposed  at  essentially 
the  same  epoch,  the  stars  will  have  virtually  identical  equatorial  coordinates. 

All  the  observations  of  a stellar  position  are  used  in  the  overlap,  and  they  are  all 
simultaneously  reduced.  Another  advantage  to  the  overlap  is  that  all  of  the  stars  act 
as  reference  material.  The  equatorial  coordinates  of  a star  are  tied  onto  the  celestial 
coordinate  system  by  the  use  of  reference  stars.  The  overlap  method  uses  the  multiple 
images  of  the  non-reference  stars  to  ‘lock’  plates  together;  while  the  reference  stars  tie 
the  whole  plate  system  to  the  sky.  This  reduction  technique,  if  carefully  and  competently 
applied,  provides  very  precise,  internally  consistent,  stellar  positions  in  the  equatorial 
frame  from  photographic  plates. 

We  then  derive  the  stellar  proper  motions  from  a comparison  of  a star’s  equatorial 
position  at  different  epochs.  In  this  study  we  have  positions  from  three  epochs,  1900, 
1956  and  1992,  allowing  a weighted  least  squares  determination  of  the  stellar  proper 
motions. 


CHAPTER  2 

MATHEMATICAL  FOUNDATION 


To  understand  the  overlapping  plate  technique  it  is  necessary  to  outline  the  theory 
of  least  squares.  The  overlap  is  essentially  a least  squares  reduction  of  stellar  equatorial 
coordinates  from  measured  rectangular  coordinates  of  those  stars  on  photographic  plates. 
This  chapter  lays  the  mathematical  foundation  for  the  development  of  this  technique. 

Traditional  Least  Squares 

Consider  a set  of  observations,  Xo,  which  we  will  assume  are  unbiased,  i.e.  without 
systematic  errors.  The  set  of  their  true  values,  x,  will  be  the  sum  of  Xq  and  a vector  e of 
corrections.  Following  our  assumption  of  unbiased  estimates,  and  assuming  a Gaussian 
distribution,  then  we  can  express  the  distribution  function  of  the  errors  as 

ip(s)  = Cexp  (2-1) 

where,  writing  all  quantities  in  vector  form,  then  formally 

a;  = i’o  + £,  (2.2) 

and  a is  the  covariance  matrix  of  Xo  (or  obviously  e). 

Assume  further  that  the  measurable  quantities,  x and  certain  parameters  a,  are 
related  by  p vector  relationships 

fx(a,x)  = 0,  (X=  (2.3) 

in  vector  form:  Fpy^i{a,x)  = 0.  The  equations  f\{a,x)  = 0 are  the  condition  equations 
(sometimes  called  “observation  equations”)  of  the  problem. 
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One  can  never  find  the  true  values  of  the  observed  quantities,  x,  but  using  the 
principle  of  maximum  likelihood,  we  can  find  an  estimate  of  x that  is  better  than  the 
original  Xo.  In  this  theory  we  seek  to  maximize  the  value  of  the  likelihood  function. 
This  occurs  when  e^a~^e  = minimum.  The  relationship  (2.3)  between  the  parameters, 
a,  and  the  target  quantities,  x,  can  thus  be  rewritten  in  terms  of  the  errors,  e,  and 
observations  Xo  as 


and  these  must  be  strictly  satisfied  at  the  solution. 

Traditional  least  squares  assumes  that  the  observations  are  of  the  same  precision, 
uncorrelated  and  that  each  equation  of  condition  contains  only  one  component  of  Xq. 
Therefore  if  the  variance  of  any  one  measurement  is  ctoo,  the  covariance  matrix  will  be 
of  the  form,  a=(Joo  I,  where  I is  the  identity  matrix  and 


This  is  the  case  in  traditional  least  squares,  where  the  values  of  e^a~^e  and 
coincide.  After  linearization  one  can  write  the  condition  equations  in  the  form 


^px  l(^,^o  “t"  — t), 


(2.4) 


(2.5) 


m 


n 


(2.6) 


written  in  matrix  / vector  form  as  Aa  + Xq  = —e.  The  parameters,  a,  occur  linearly 
in  the  condition  equations,  and  each  equation  contains  exactly  one  observation.  (In 
general  any  number  of  observations  may  occur  in  any  of  the  condition  equations.)  It  is 
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easy  to  find  the  normal  equations  that  give  estimates  of  a 

Aa  + Xo  = —e 

(2.7) 

— ► a^A'^Aa  + 2a^A^x<>  + xjxo  = e^e. 

(Note:  e^e  and  a"*'A^jCo  are  numbers  and  we  may  transpose  any  of  their  terms  without 
changing  their  values.)  Now  differentiate  with  respect  to  each  component  of  a and  set 

it  equal  to  zero  to  find  the  maximum: 

d(e^e)  d{aT  A'^  Aa + 2aT  A'^Xo x^ x„) 

da  da  ^2.g) 

= 2(A'^Aa  + A'^Xo). 

If  the  right  hand  side  of  (2.8)  is  set  equal  to  zero,  then  the  condition  becomes 
A^^Aa  = — A^Xo.  By  rearranging,  we  get  the  traditional  least  squares  solution  in 
its  simplest  form 


a = — (A^A)  ^A^Xo- 


(2.9) 


Generalization  of  the  Traditional  Method 

We  can  use  a more  general  least  squares  treatment  by  dropping  some  restrictions. 
Following  the  procedure  of  Brown  (1955)  (see  also  Jefferys,  1980,  1981;  Eichhom, 
1993),  we  will  make  use  of  Laplace  multipliers,  A.  At  the  solution  the  condition 
equations,  F,  equal  zero.  Therefore  the  value  of  S = e^a~^e  is  the  same  as  the  value 
of 

s*  = - 2A'^F(a,  x<,  + c).  (2.10) 

Therefore  the  values  of  S and  S*  reach  their  minima  at  the  same  values  of  e and  a. 
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The  parameter  vector,  a,  has  n components  and  e has  m components,  so  these  two 
give  (m+n)  free  components.  The  vector  of  condition  equations,  F=0,  has  p components. 
Since  F is  restricted  to  = 0,  there  are  only  (m+n-p)  free  components.  Differentiate  S* 
w.r.t.  the  components  of  e and  a,  thus 

j =2<r-^e-2X'^A 


(2.11) 


-2A'^A 


where  we  have  used  the  abbreviations 


and 


• 

U • 

. Ml.} 

dt„ 

f =Xo,a  = (3o 

/i/j.  . 

' dai 

da„  ' 

' doi 

. Mjlj 

dan  ' 

(2.12) 


(2.13) 


At  the  solution,  the  following  equations  must  thus  be  rigorously  satisfied; 

a~^e  = X^A  m equations 

A^yl  = 0 n equations  (2.14) 

F = 0 p equations. 

The  vector  A has  p unknowns;  from  (2.14)  there  are  exactly  as  many  unknowns  as 
there  are  equations. 

These  equations  are  rigorous  at  the  solution  so  A and  X must  be  evaluated  at 
the  solution  of  a and  x.  These  equations  are  nonlinear  and  solved  by  an  iterative 
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procedure  (Jefferys,  1981).  Assuming  some  approximate  value,  Gq,  for  a and  using  as 
approximations  to  x the  observations,  ;Co,  we  write  a=ao+a  and  expand  F as  a Taylor 
series. 

Formally, 

a = a„  + a,  x = z<,  + £ 

F(x,  a)  = F(xo,  ao)  + X„£  + A„a  + 0[2]  (2.15) 


where  Xo  and  Aq  from  equations  (2.12)  and  (2.13)  and  evaluated  at  Xo  and  ao.  From 
the  first  of  equations  (2.10)  we  get  e = crX^A;  insert  the  expansion  of  F(x,a)=0: 


F = 0 = Fo  + Xof  + AqQ: 


— ► Fq  F XotrX*^A  -j-  AoQ^  — 0 


(2.16) 


where  Fq=F{  XqAo)  combine  with  the  second  of  equations  (2.10);  A^A  = 0 so  that 
we  may  write 


fX<rX'^ 

V A-r 


/o\ 

U/ 

(2.17) 


These  are  the  generalized  normal  equations.  Note  that  we  need  A to  find  e,  but  if 
X^crX  is  nonsingular  then  we  can  eliminate  A and  find  the  corrections  to  the  parameters, 
a.  This  may  be  done  as  follows: 

F<,  + XcrX^A  + Aq  = 0 

- (XcrX'^)“ Vo  + A + (XaX'^)' Vq  = 0.  (2.18) 


Hence 

A = -(XaX'^)"  Vq  - (X<rX'^)"Vo 


^ A = -W(Aq  + Fo) 


(2.19) 
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where  W = (XcrX"^)  but  since  A"*- A = 0 we  get 


A'^(X(tX'^)  Va  = -A'^(X(tX'^)  Vo 


a = — 


-1 


A’^(X<tX'^)  V ‘a'^(X(tX'^)  Vo 


(2.20) 


This  is  a more  general  solution  than  that  of  equation  (2.9). 
The  formal  errors  e of  the  observations  are  then  given  by 


e = (tX'^A  = -aX'^W(Aa  + Fo).  (2.21) 

It  can  be  shown  (cf.  Brown,  1955;  Jefferys,  1980,  1981)  that  (A^WA)  is  the 
covariance  matrix  of  a,  so  that  the  square  roots  of  its  diagonal  terms  are  the  standard 
deviations  of  the  corresponding  components  of  a. 

It  is  instructive  to  compare  the  results  of  the  traditional  approach  with  that  of  the 
more  general  approach.  In  the  simplified  case,  where  we  use  noncorrelated  observations 
of  equal  precision  and  when  exactly  one  observation  occurs  in  each  equation,  X=-I. 
The  corrections  to  the  parameters  simplify  to  a = -[A’^A]~^A’^Fo.  The  observation 
errors  also  simplify  to  e = |l  - A[a'^A]”^A'^|Fo.  both  of  which  are  traditional 
least  squares  solutions,  independent  of  the  variances.  Also  from  comparison  we  can 
see  that  the  matrix  (XcrX"^)  ^ has  the  same  effect  as  a weight  matrix  in  the  traditional 
approach. 

This  matrix,  (X<rX^)  is  not  always  nonsingular,  as  for  example  if  there  are 
equations  in  F that  have  no  observations  but  act  as  parameter  constraints.  This  will 
not  occur  in  our  problem. 
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In  the  general  case  the  usual  scheme  is  to  set  up  iterations  as  follows: 


fX.aXj  A,\ 

V Aj  0 ^ 

^ A^+1  \ / F^'\ 

+ 

W+i)  \0  J 

/o\ 

\o) 

where  X,  = ^ ^ ^ A = 

VlaJ  a = Co  + cx,/  ' 

x = Xo  + £ X = Xo  + e 


(2.22) 


This  does  not  always  converge,  and  more  sophisticated  schemes  are  usually  needed 
to  ensure  convergence.  For  this  purpose  the  important  point  is  that  in  this  scheme  the 
iterations  are  performed  with  both  the  e and  the  a as  variables.  We  will  apply  this 
theory  to  derive  the  parameters  that  define  the  position  of  a stellar  image  on  a single 
photographic  plate. 


CHAPTER  3 

FINDING  STELLAR  COORDINATES 


Now  apply  the  theory  of  least  squares  to  the  determination  of  equatorial  coordinates 
from  the  rectangular  coordinates  of  a stellar  image  on  a photographic  plate.  The  first 
section  presents  the  reduction  for  single  plates,  following  Eichhom  (1984). 


Consider  the  data  from  a single  photographic  plate.  On  it,  there  will  be  m 
reference  stars  and  n stars.  The  position  estimates  for  the  reference  stars  are 

available  from  some  catalogue,  with  their  respective  variances,  pi^co^Spc  and  cti/.  First 
make  the  (reasonably  correct)  assumption  that  the  aye  and  Syc  are  uncorrelated.  The 
Astrographic  Catalogue  (or  direct  measurements  from  a photographic  plate)  provides 
measured  rectangular  coordinate  estimates  Xy  and  yy.  The  variances  v and  <f>  for  these 
coordinate  estimates  must  be  obtained  somehow. 

Let  us  introduce  “standard  coordinates,”  ^ and  tj,  to  facilitate  the  reduction.  These 
are  related  to  the  right  ascension  and  declination.  The  relationship  between  ^ and  tj 
on  one  hand,  and  a and  S on  the  other  hand  describes  the  geometry  of  a gnomonic 
projection: 


Positions  from  Single  Plates 


(3.1) 
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or  inverted 


cos  (5  sin  (q  — ao) 

sin  6 sin  6o  + cos  6 cos  So  cos  (a  — ao) 


(3.2) 


sin  6 cos  ^0  ~ cos  6 sin  cos  (a  — ao) 
sin  6 sin  ^o  + cos  S cos  Sq  cos  (a  — ao) 


Consider  the  ideal  situation  where  the  optics  of  the  telescope  provides  exactly  the 
same  imaging  as  a pinhole  camera.  The  rectangular  coordinates  Xi/.yj/  would  thus  be 
related  to  these  standard  coordinates  purely  by  the  focal  length  and  the  orientation  of 
the  plate  in  the  focal  plane.  However,  a telescope  never  provides  an  ideal  gnomonic  (i.e., 
pinhole  like)  projection.  The  deviations  of  the  telescope’s  properties  from  producing  an 
ideal  gnomonic  projection  — these  are  called  aberrations  — produce  a variety  of  optical 
effects  that  affect  the  position  of  a star  image.  Even  if  the  telescope  did  produce  images 
exactly  in  the  manner  of  a pinhole  camera,  we  would  still  rely  on  experimental  estimates 
of  the  telescope’s  parameters.  For  example  we  need  the  focal  length,  the  location  of 
the  plate  center  and  the  orientation  of  the  plate  with  respect  to  the  optical  axis.  All 
the  deviations  and  telescope  parameters  are  included  in  the  model  and  estimated  from 
the  observations. 

In  practice,  we  assume  the  relationship  between  the  standard  coordinates  and  the 
rectangular  coordinates  of  the  i/-th  star  to  be  of  the  form 


f Xu  '' 

— s 

\yu) 

U-/ 

where  s is  an  accurate  approximation  of  the  focal  length  of  the  telescope,  is  the 
model  matrix  (with  two  rows),  and  a is  a vector  of  parameters  (the  parameter  vector). 
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The  values  of  the  parameter  vector  a can  be  found  from  a linearized  least  squares 
analysis  if  |S^a|  <<  |s  1 [. 


The  model  matrix,  S,  and  the  concomitant  parameter  vector,  a,  can  take  many 
forms  depending  on  the  properties  of  the  telescope.  For  this  example  we  will  use  the 
simple  “six-constant  model.”  In  this  model  relationship  (3.3)  takes  the  form 


X 

y 


n T)  \ Q Q Q\ 

Vo  0 0 ^ 7?  ly 


6 

c 

a 


(3.4) 


b' 

\c) 

where  each  of  the  parameters  a ...  c has  a physical  significance. 

Parameters  c and  c'  account  for  differences  between  the  origins  of  the  measurable 
rectangular  image  coordinates  (x,y)  and  the  standard  coordinates  (^,7]).  a and  b'  correct 
the  focal  length.  Rotation  of  the  plate  is  corrected  for  by  the  parameters  b and  a'.  This 
relationship  (3.4)  models  any  effect  which  produces  a linear  relationship  between  (x,y) 
and  systems. 

For  this  analysis  we  use  a matrix  that  contains  up  to  16  independent  parameters, 
but  the  underlying  reduction  remains  in  principle  the  same.  Obviously,  larger  and  more 
sophisticated  models  may  account  for,  and  correct  more,  complicated  errors. 
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Use  of  the  relationship  (3.3)  will  render  the  condition  equations  for  the  reference 
stars  to  be  of  the  form 


Note  that  the  are  matrices  of  dimension  (2xm),  where  m is  the  number  of  components 
of  vector  a. 

In  equation  (3.5),  the  xi,yi,...,Xm,ym  are  the  observations,  and  the  ai,<5i...,Qm,<5m  and 
the  components  of  the  vector  a are  the  unknowns.  For  linearizing  these  equations,  we 
need  approximate  values  for  the  Qi,...,<5m.  An  analysis  of  a single  plate  only  incorporates 
reference  stars  (i.e.,  stars  for  which  estimated  a and  6 are  available).  The  catalogue 
position  estimates  of  these  reference  stars  provide  the  additional  condition  equations 
which  are  necessary  to  make  the  system  determinate: 


= 0. 


(3.5) 


^ \ / (<^m , ^ 


\ \ 2/m  / V 7m  ( ® m 1 ) / 


/ (ttlc  - «l)cos^i 


G = ^2c  — ^2 


= 0. 


(3.6) 


i^mc 


\ ^mc 
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Here,  qj/c— ,^i/c  are  observations  (coordinate  estimates  from  a reference  catalogue), 
while  au,...,Su  are  the  same  unknowns  as  in  H=Q  from  (3.5). 

We  combine  all  the  condition  equations  into  one  set: 


F = 


(3.7) 


F has  the  same  function  as  in  the  previous  analysis  of  the  least  squares  solution.  From 
equation  (2.20)  the  corrections  a to  the  parameters  will  be  given  by 


A'^(X(tX'^)  ^a]~\'^(X<tX'^)'V„ 


(3.8) 


where  the  i/-th  element  of  each  matrix  is  X = (4^)  , , , A 

V dx  /a=ao+Q'i/,J:=io+e 


(dF\ 

V da  / a=a. 


o+au,x=Xo+£ 


, Fo  = F{Xo,ao). 


All  the  equations  contain  exactly  one  observation;  therefore,  we  have  the  (classical) 
happy  circumstance  that  = (^)  = I,  the  identity  matrix.  Our  first  assumption 
that  the  observations  and  catalogue  coordinates  are  uncorrelated  gives  us  a covariance 
matrix  of  the  form 

i > 

0 


Fm 


Pi 


(3.9) 


<^i 


Pm 


\ 


Om) 


let  cr  = 


(Xxx  0 


0 ^ act 

Equation  (3.8)  then  simplifies  to 


a = -[A’^<r  ^A]  ^A^(x~^Fo. 


(3.10) 
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For  an  initial  approximation  we  take  the  components  of  the  parameter  vector,  a,  to  be 
zero.  Now  insert  in  the  left  hand  side  of  the  condition  equations  (3.5)  the  observed  pairs 
x,y,  put  a=0  and  calculate  the  ^ and  tj  from  the  catalogued  a and  S of  the  reference  stars. 
The  right  side  will  not  be  equal  to  zero  but  a vector  d with,  we  hope,  relatively  small 
components:  the  residual  vector.  Also  let  the  corrections  to  the  positions  be  a vector  fi 
defined  by  small  position  increments  (cosSda,  d6).  These  are  one  set  of  unknowns. 

Equation  (3.8)  thus  formally  becomes  F = = (O’  ^ residual 


vector. 


With  O'approx — ^catalogue  ^tid  ^approx — ^catalogue  ^he  residual  vector  for  G will  be  zero 


by  definition.  The  corrections,  q,  to  the  parameters  become  a = 
vector  of  corrections  to  the  reference  star  positions,  i.e. 


, where  /3  is  the 


/ cos6idai  \ 
dbi 

cos  dctifYi 

dSffi 

V 

and  a is  the  vector  of  the  plate  parameters. 
Equation  (3.8)  now  simplifies  to 


(d\ 

= — [A^(t  ^A]  ^A’^<7~^ 

\a  y 

\o) 

(d\ 

[A'^^-^A] 

\aj 

koJ 

(3.11) 


(3.12) 
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First  we  must  investigate  the  structure  of  A. 
We  have  defined  A as 


A = 


f^\  ^ (d{H,G) 
\da  J V o) 


Consider  the  terms  in  equation  (3.12) 


V 0 ) 


Multiply  this  by  the  matrix  A: 


A^a  = 


+ (Tea  sB'^o-^iS ' 


Substituting  this  into  equation  (3.8),  we  get 


s^B'^orJj?  B + sB'^CTrr  2 ' 
B 


S't.t-'S 

“ ^ XT  ^ 


.aJ 


( sB'^cr,.,.*  (T 


XX  ^ XX 


“ XX 


(d\ 


\oy 


(3.13) 


(3.14) 


(3.15) 


(3.16) 
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expanding 


+ o-oa)/?  + sB'^aJZa  = sB'^<7~^d 


Solve  for 


5STariB^+  ET<TriSa  = 


/?  = (s'B'^cTj.^^B  4-  cTao)  ^sB’^cr„^(d  - Sa) 


and  a 

sE^a~^B(^{s^B'^cr~^B  + ffaa)~^ sB'^cr~^(d  - 'Eafj  + S'^aJ^Sa 

— wT.-lj 

- “ «^xx  « 

whence 

“^{<^xx  - «*<^jiB(s^B'^«rJ^^B  + ct„o)"‘b'^(t-^‘ |sa  = 

- «^<^jiB(s^B'^<T"jB  + (Tc.c.y^B'^ar^^'^d, 

SO  that, 

S^JSa  = S^Jc/  or  a = (E^JH)  ^S'^Jd 

where  J = a~^  - + craa)  ^B'^a-}. 

This  matrix  J can  be  simplified  using  the  matrix  inversion  lemma 

J = (cT,,,  + Ba;iB'^)"' 

which  is  a block  diagonal  consisting  of  (2  x 2)  blocks  on  the  main  diagonal, 
parameter  vector,  a,  more  explicitly  is 

m 


(3.17) 

(3.18) 

(3.19) 

(3.20) 

(3.21) 

(3.22) 
The 


a = 


(3.23) 
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Now  we  must  find  more  explicit  expressions  for  the  matrix  J.  We  could  go  directly 


formula  is  simple  but  requires  first  approximations  for  the  spherical  coordinates,  a, 6.  It 
is  better  to  express  the  matrix  B in  terms  of  quantities  we  already  have  estimates  for, 
e.g.,  and  rj  for  which  first  estimates  are  just  -jx  and  ^y.  In  the  case  of  the  astrographic 
catalogue,  estimates  of  the  parameters  give  very  close  estimates  for  the  ^ and  ri  in  terms 
of  X and  y. 

Eichhom  (1985)  showed  if  we  define  the  following  quantities. 


to  equation  (3.22)  but  this  would  entail  evaluating  B.  To  evaluate  B from  the  initial 


5=  {^R-sm6o)/T 


(3.24) 


U = -f  1)  cos  6o))/T 


V = (i2‘(cos^o  — r)s\nSo))/T 


W = {{(^Rcos6o)IT)  + VlR 


we  may  write  the  matrix  B = ^ ^ ^ f ^ . If  one  further  introduces  the  quantities 

Q = -UWp^SVa 


Y = U^p  + l/V 


(3.25) 


Z=  W'^pYS'^a 

where  p = cr a and  a = as . Then  from  equation  (3.22)  , 


(3.26) 
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Finally  we  can  also  use  the  quantities  above  to  simplify  the  corrections  to  the 
positions,  again  from  Eichhorn  (1985), 


^ pa  + Y u + Z <p) 


(3.27) 


Example 

This  completes,  in  essence,  the  theoretical  analysis  of  a single  plate.  It  is  helpful  to 
summarize  what  the  steps  are  in  extracting  the  position  estimates  from  an  astrographic 
plate  as  in  the  immediate  problem. 

Combine  the  published  spherical  coordinates  of  the  reference  stars  (qc,Sc)  w'ith 
the  published  plate  center  (oo.^o)  in  the  trigonometric  relations  (3.2)  and  calculate  the 
standard  coordinates  We  use  the  estimated  plate  parameters  from  the  AC  to  find 

approximate  rectangular  coordinates  (x',y')  for  the  equatorial  coordinates  of  the  reference 
stars.  Compare  the  measured  coordinates  (x,y)  of  the  AC  and,  using  the  magnitude  as 
an  independent  check,  identify  the  reference  stars. 

After  this  process  we  have  measured  coordinates  (jc,y)  matched  with  the  reference 
stars’  spherical  coordinates,  (oc^c).  Also  for  each  plate  we  have  initial  estimates  of  the 
center  and  the  parameters.  The  next  step  is  to  find  approximate  standard  coordinates 
using 


where  the  parameters  represent  any  shift  of  the  zero  point  from  the  initial  estimates.  We 
use  these  approximate  values  (^',1]')  for  no  other  purpose  than  to  calculate  the  matrix 
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J from  equation  (3.26)  and  define  the  matrix  2 in  equation  (3.23).  The  approximate 
values  (^',7/0  are  fully  sufficient  for  these  calculations. 

Standard  coordinate  estimates  are  found  from  equation  (3.2)  using  the  spherical 
coordinates  of  the  reference  stars  with  the  assumed  plate  centers.  We  use  these  values  to 
calculate  the  residual  vector,  d.  If  we  assume  an  initial  value  of  zero  for  the  parameters, 
then 


Xi  -S.^1 

yi  - s.T]i 


d = 


(3.29) 


. Vm  . 


and  from  equation  (3.23),  for  the  reference  stars  on  the  plate. 


a = 


(3.30) 


we  find  the  vector  a.  As  we  have  assumed  initial  approximations  in  the  residual  vector 
to  be  zero,  this  a will  be  the  actual  parameters,  except  a(3)  and  a(6)  will  be  corrections 
to  the  originally  assumed  values. 

We  have  found  the  single  plate  parameter  estimates.  Using  these  we  can  estimate 
from  the  rectangular  coordinates  (x,y)  published  in  the  AC,  the  standard  coordinates 
of  the  field  stars.  These  in  turn  provide  spherical  coordinates  using  the  inverse  of  the 
gnomonic  relation.  We  can  now  use  the  overlapping  plate  technique  to  improve  on 
these  first  estimates  of  the  field  stars’  equatorial  coordinates. 
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Positions  from  Overlapping  Plates 
Solution  with  Unconstrained  Parameters 

Single  plate  solutions  do  not  use  all  the  available  information.  The  most  important 
information  not  used  is  the  constraint  that  the  same  star  on  two  plates  of  the  same  epoch 
must  have  the  same  set  of  equatorial  coordinates,  because  a star  can  only  be  at  one  place 
at  one  time.  In  a single  plate  solution,  the  data  on  each  plate  give  independent  “best” 
estimates  that  can  be  then  averaged  to  give  a “combined  best”  (even  “bester”?)  estimate. 

This  process  leads  to  systematic  errors  because  the  available  plate  parameter  esti- 
mates are  not  the  true  ones  (Eichhom  and  Wiiliams,  1963).  The  individual  parameter 
errors  generate  systematic  errors  in  the  star  positions  calculated  with  them.  This  is 
aggravated  by  ignoring  the  above  constraint;  enforcing  more  constraints  generates  more 
accurate  parameters.  An  overlap  solution  enforces  this  constraint,  producing  more  ac- 
curate plate  parameters,  thus  inevitably  leading  to  smaller  systematic  errors  in  the  star 
positions.  The  experience  of  many  investigators  has  shown  however  that  the  reduction 
model  (i.e.  the  matrix  S)  must  be  carefully  chosen  to  be  the  most  realistic  available. 
Assume  that  there  are  altogether 

1.  n at  least  partially  overlapping  plates,  which  are  numbered  i/=l,...,n, 

2.  m stars  in  the  region  which  are  either  on  more  than  one  plate  or  a reference  star; 
H is  the  current  star  number, 

3.  nix  of  the  m stars  in  the  region  are  reference  stars, 

4.  The  numbers  of  the  reference  stars  are  I-Iti,  ■■■■, 

5.  mj/  stars  on  the  ui\\  plate. 
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The  measures  (estimated)  rectangular  coordinates  of  the  /ith  star  on  the  i/th  plate  are 


y/zj/- 


The  equations  of  conditions  are  identical  with  those  in  a single  plate  solution.  The 
first  step  in  the  overlap  is  to  find  estimates  of  the  spherical  coordinates  a and  8 of  the 
field  stars  from  single  plate  reductions  and  average  those  that  appear  on  more  than  one 
plate.  This  yields  the  required  estimated  equatorial  coordinates  for  all  stars  that  are 
in  turn  used  to  calculate  the  standard  coordinates.  There  must  be  no  loss  of  significant 
figures  in  this  process. 

It  is  essential,  at  this  point,  to  use  the  same  estimates,  a^ofipo,  for  calculating  the 
standard  coordinates,  of  each  star.  For  the  reference  stars,  the  catalogued 

estimates,  acfxficp,  should  be  used  as  initial  approximations  to  these  stars’  equatorial 
coordinates. 

The  condition  equations,  generated  by  the  measurable  coordinates  x,y  for  the  i/th 
plate,  will  therefore  be  of  the  form 


/ 


H = 


\ ynt,^v  — j 


V 


/ 


= 0 


(3.31) 


Note  that  equations  (3.31)  will  be  satisfied  only  for  the  true  values  of  the  measurable 


coordinates  x and  y and  the  final  estimates  for  the  plate  parameters  Ui/. 
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In  these  equations,  we  have  assumed  that  the  (most  likely  nonconsecutive)  numbers 
of  the  stars  on  the  z/th  plate  run  from  Hui,  ■■■,  which  necessitates  the  complicated 
indexing.  Note  that  each  /xj/  is  just  one  natural  number;  an  image  of  the  same  star  on 
plates  1/  and  i>'  will  yield  two  numbers  k and  A such  that 

Assign  the  numbers  /x  to  the  stars  in  some  organized  fashion,  for  example  in  order 
of  increasing  right  ascension.  The  stars  next  to  each  other  are  then  likely  to  be  measured 
on  plates  that  overlap  one  another  at  least  partially.  This  will  be  elaborated  on  in  the 
example  given  below. 

Concerning  the  reference  star  condition  equations,  it  is  best  to  write  the  equations 
corresponding  to  G=0  for  all  reference  stars  together.  Otherwise  if  we  write  G after  the 
equations  (3.31),  it  will  require  an  awkward  criterion  to  make  sure  than  the  equation 
Gfi=0  is  set  up  only  once  for  each  reference  star. 

The  measured  coordinates  of  the  stars’  images  on  each  plate  provide  the  set  of 
equations 


1/  = 


/H,\ 

th 

\h„J 


= 0. 


(3.32) 


The  reference  stars,  that  is  those  that  have  estimates  Oc^c  of  their  spherical  coordinates 
listed  in  some  catalogue,  generate,  in  addition  to  H=0,  the  following  set  of  condition 
equations 

^CUrl  ~ + f/irl 

~ + ^Mr3 


G = 


(Qcu  —Oil  ) COS^I,  + £„ 

\ ^rn\r  f t^rrtir  ' ^ t^rlrvur 

\ ^C^rmr  ^t^rmr  "h  / 


= 0 


(3.33) 
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where  c is  a catalogue  value.  These  equations  reflect  the  previous  statement  that  the 
numbers  of  the  reference  stars  are  /^n,  ,/^r/wr  • 

There  is  an  essential  difference  between  equations  (3.31)  and  (3.33).  Equation 
(3.31)  holds  for  the  “true”  values  of  the  observables,  i.e.  the  rectangular  coordinates  of 
the  stellar  images  on  the  plates.  The  equations  are  satisfied  rigorously  after  we  have 
found  the  “definitive”  values  (estimates)  for  the  plate  parameters  aiy. 

Equation  (3.33)  is  somewhat  different.  Instead  of  observables  (the  spherical  co- 
ordinates themselves),  we  have  used  their  catalogued  estimates;  the  observables  are, 
therefore,  also  adjustment  parameters.  To  be  able  to  equate  (3.33)  to  zero  we  need  to 
add  corrections  e and  ^ to  the  appropriate  catalogue  estimates.  Normally  the  numbers 
/^rk./irk+i  of  two  “neighboring”  reference  stars  in  equations  (3.33)  will  not  be  consec- 
utive. In  typical  situations  only  a few  of  the  stars  on  each  plate  will  be  reference 
stars,  and  they  will  usually  not  have  consecutively  ordinal  numbers.  Again  this  will 
be  illustrated  by  the  example. 

Equations  (3.31)  and  (3.33)  combine  to  form  the  condition  equations; 
F = ^ ^ ^ =0.  Since  only  one  observation  occurs  in  each  equation  the  matrix  X 
(being  the  Jacobian  matrix  of  the  condition  equations  F = 0,  with  respect  to  the  ob- 
servables, x)  will  be  the  identity  matrix.  This  simplifies  our  least  squares  solution, 
such  that. 


= I 


(3.34) 


therefore. 


(3.35) 
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From  equation  (2.20),  the  correction  to  the  parameters,  a,  is  given  by  a = 
(XcrX"^)  A'^(X(7X’^)  ^ Fo.  When  X = I this  simplifies  to,  a = 


The  next  step  is  to  find  the  matrix  A = (^),  where  a are  the  adjustment 
parameters.  Let  ^ be  the  corrections  to  the  stars’  spherical  coordinates  and  a the  plate 
parameters,  thus  a = example  |o|  <<  1 and  if  the  parameter  vector  is 

a = 0 then,  the  values  we  obtain  from  equation  (2.20),  will  be  the  actual  parameter 
values  rather  than  the  corrections. 

Since  A is  the  Jacobian  matrix  of  the  condition  equations  with  respect  to  the 
adjustment  parameters,  we  have 


A = 


\d{(3,a)  ) 


From  equation  (3.32)  we  have 


dH, 

da 


/M^\ 

da 

an, 

da 

\ da  ' 


and  for  any  Hi/, 


(3.36) 


(3.37) 


dH,  ^ dHA 

da  \ dd  ’ da  ) 


(3.38) 


In  this  equation  ^ is  the  vector  of  star  parameters,  {a^  cos  6^,  ^^),  and  a is  the  vector 
of  plate  parameters,  taken  one  plate  after  another:  = (af,  aT,  We  consider 

the  condition  equations  Hi/=0  which  arise  from  the  measured  rectangular  coordinates 
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of  the  star  numbers,  on  the  vih  plate.  Break  down  Hj/  further  and  write 


H,  = 


(H,\ 

Hz 

\hJ 


(3.39) 


where  is  a vector  of  dimension  (2  x 1). 

d H 

Since  there  are  altogether  m stars,  the  Jacobian  matrix  —5^  is  of  dimension  (2 
X 2m):  = (0 . . . . . . O)  where  the  0 are  matrices  of  dimension  (2  x 2) 

and  B„^  p = ( -rr — — r V are  also  of  dimension  (2  x 2).  The  reason  for 
the  particular  structure  of  j is  that  each  pair  ^,7/  depends  on  only  one  pair  a,^. 

We  get  therefore  from  “each  plate” 


= Bp.  (3.40) 


0 . 

0 

0 . 

0 

...  ox 

0 . 

0 

0 . 

• Bu  1/ 

0 . 

0 

...  0 

dHp 

-s 

. 

. 

dp 

• 

• 

• 

• 

• 

• 

• 

\0  . 

0 

0 . 

0 

0 . 

...  0/ 

This  matrix  is  of  dimension  (2m//  x 2m).  Each  line  contains  exactly  one  (2  x 2) 
matrix  that  is  not  a null  matrix;  and  each  column  at  most  one  (2  x 2)  matrix  that  is  not 
a null  matrix,  namely  when  the  an  image  of  the  star  is  measured  on  the  uth  plate.  One 
could  say  that  is  then,  and  only  then,  not  a null  matrix  when  /i  e i-C-,  when 


the  number  //  of  the  star  is  among  those  whose  images  were  measured  on  the  uth  plate. 

/Bi\ 

then  we  may  wnte 


If  we  introduce  B = 


\Bn/ 


dH 

dp 


= -sB. 


(3.41) 


The  matrix  B obviously  has  the  dimensions  (2  ^ x 2m.) 

i/=i 
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/ V \ 


If  we  write  Sj,= 


•‘HviV 


, which  has  the  dimensions  {2mu  x ly)  where  ly  is 


the  number  of  components  of  the  vector  Uy,  then  we  see  from  equation  (3.31)  that 


^ = (0  ...  0 -s, 

dH  _ dH^ 
da  d{ai,a3,  ■ • • , a„) 


0 ...  0)  and 


dH, 

dH. 

dH,. 

dat 

doi 

' ‘ dOi 

dH, 

dco 

dHo 

da. 

dH,. 
' * da^ 

dH, 

dH,. 

dH,. 

day 

' * day 

/2.  0 

0 S, 

Vo  0 


0 \ 
0 


2„/ 


(3.42) 


We  further  see  that 


((q>  — a)  cos  6 ^ 

Sc-S  / 

d/3  d(a  cos  S,  6) 


(3.43) 


where  K is  a matrix  of  dimension  (2nir  x 2m).  It  has  exactly  one  I2  in  each  double  line 


and  at  most  one  I2  per  double  column,  namely  for  those  star  numbers  that  correspond 


to  reference  stars.  We  see  therefore  that  the  structure  of  K is 

/O  ...  OI2O  ...  000  ...  000  ...  0\ 

0 ...  000  ...  OI2O  ...  000  ...  0 

K = . . . . 

• • • • 

• • • • 

\0  ...  000  ...  000  ...  OI2O  ...  0/ 


(3.44) 


Since  G is  independent  of  the  plate  parameters  we  see  that  = 0 and  of  dimension 

n 

(2mr  X ^ p„). 

The  matrix  A contains  the  same  derivatives  in  an  overlap  solution  as  a single  plate 
reduction,  but  its  structure  is  more  complicated  for  an  overlap  reduction.  Generally,  A 
can  be  written  in  the  form 


-A  = 


sB  S 
K 0 


(3.45) 
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which  is  analogous  to  its  form  for  a single  plate  reduction.  The  deeper  structure  of  A for 
an  overlap  reduction  deserves  close  examination,  and  will  be  analyzed  in  the  example. 

The  matrix  B will  be  made  up  of  individual  blocks  of  dimension  (2  x 2),  B^j/,  at 
each  position  where  the  /ith  star  appears  on  the  i/th  plate.  Each  plate  will  generate  my 

n 

separate  lines;  size  2 ^ x 2m. 


v=\ 


Consider  the  normal  equations  (A^<j  ^A)a  = — A^cr  that  is  to  say 


a 


-A^a-^Fo. 


(3.46) 


Let  us  therefore  investigate  the  structure  of  the  term  (A^cr  ^A).  Multiplying  out 


A’^ct  ^A  = 


E + ...  sBTa-iS„\ 

5SJ'c7"^Bi  ...  0 


...  S^a-^Sn/  (3.47) 


A^cr  ^A 


L sB^a-^S 


with  L = + K'^ct-^K. 

Consider  the  term  (B^cr“^B),  and  the  discussion  above  of  the  structure  of  B;  on 
each  row  there  is  one  line  pair  for  each  star,  consisting  of  a (2  x 2)  block,  otherwise 
that  row  will  contain  only  zeros.  On  each  column  there  will  again  be  a (2  x 2)  block 
corresponding  to  a star  that  appears  on  the  plate  corresponding  to  the  row  it  appears  in. 
(The  diagonal  covariance  matrix  a will  not  affect  the  structure  of  B^B.) 

The  matrix  B^B  is  a multiplication  of  two  matrices  of  dimension;  (2m  x 2 E) 

V=1 

n 

and  (2  x 2m),  resulting  in  a (2m.  x 2m)  matrix.  The  /<th  row  and  column  pair 

i/=i 

n 

will  have  an  (2  x 2)  block  equal  to  E) 

j/=i 
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Consider  the  product  sB^cr~^2;  again,  the  covariance  matrix  will  not  affect  the 
structure.  The  individual  elements,  s'Bjcr^^Zt,,  are  sparse  matrices  of  dimensions 
(2m  X l^).  Nonzero  elements  of  the  product  will  occur  only  where  the 

/ith  star  appears  on  the  uth  plate.  Thus  the  structure  of  this  matrix  will  be  similar  to 
that  of  with  the  j/th  column  pair  replaced  by  a line  of  width  Ij/  , i.e.  is  a 
(2m  X 2m)  matrix. 

Finally,  the  product  is  block  diagonal  with  the  individual  blocks  equal 

to  the  products  of  each  plate’s  model  matrices  and  covariance  matrices.  This  ends  the 
consideration  of  the  left  hand  side  in  (3.46). 

n 

Consider  the  right  hand  side  of  (3.46).  This  will  be  a vector  with  2m  + 


components 


/ " \ 


i/=i 


wT.-l  j 
“1  “1 


(3.48) 


/ 


We  arrive  at  the  solution  by  first  eliminating  the  star  parameters  and  solving  for 
the  plate  parameters  only.  Consider  (3.46):  = -A^cr~^Fo.  We  can 

substitute  A^cr“^A  and  A^a~^Fo  from  equations  (3.47)  and  (3.48)  to  find  that 


(3.49) 


or  multiplying  out 


L(3  + sB^(t  ^2a  = sB^(r~^rf 


^B/3  + S'^cr  ^Sa  = 


(3.50) 
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The  first  row  gives 


/3  = sL-'B'^aZ^id-Ea) 
and  /?  = si  j 

\ V=1  / „=l 


(3.51) 


The  vector  /3^=  (da^cos  S^i,  is  different  from  the  quantity  ^u-  The  term 

is  a Kronecker  symbol  which  will  be  equal  to  zero  if  the  //th  star  is  not  a reference 
star  and  equal  to  one  if  it  is  a reference  star.  The  components  of  the  matrix  ^ are  the 
corrections  to  the  star  coordinates  individually  and  amount  to  weighted  means  of  all 
the  frames  on  which  a particular  star  occurred. 

Substituting  this  into  the  second  row  of  (3.50)  gives 

s^S^(T~^BL“*B^(T“^(d  — Sa)  + S^<r“^Sa  = 2’^cr“^d 


^ S'"(<t;^  - sV;^BL-B^<T;^^)(d  - 2a)  = 0 


(3.52) 


simplifying 

a=  (2T"J'2)''2Vc/ 


(3.53) 


where  J'  = a~^  — s^cr^ 

The  inversion  lemma  cannot  be  used  on  J'  because  the  matrix  kTct“^K  in  L is 
singular.  The  matrix  J'  is  not  block  diagonal  but  assumes  the  pattern  of  B^B  because  L 
is  block  diagonal.  The  product  S^J'Sa  is  ’banded- bordered’  and  there  exist  efficient 
routines  to  invert  it. 

It  is  instructive  to  examine  the  form  of  the  Laplace  multiplier  A introduced  in 
equation  (2.10).  From  equation  (2.9) 


A = -W(Aa-fFo) 


(3.54) 
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where  W = (XcrX"^)  In  the  overlap  case 

X = I 


.-.w  = {XcrX'^y^  = a~\ 

Substituting  this  into  equation  (3.54),  we  get 


A = -(T~^(Aa  + Fo)  = cr~^A 


+ (T  ^Fo, 


(3.55) 


(3.56) 


and  substituting  for  a,  A and  Fo 


From  equation  (3.5)  the  formal  errors  in  the  observations  are  therefore 


c = crX'^A  = 


/ sB/J  + Ha  — d\ 

V K/?  ;• 


(3.57) 


(3.58) 
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Example 


Plate  2 

4 

O 

3 

o 

0 

5 2 

ignored 

0 0 

Reference  Star 


Field  Star 


0 

ignored 


Plate  1 


Figure  2:  Arrangement  of  Stars 


As  in  the  single  plate,  it  is  helpful  to  summarize  what  the  steps  are  in  evaluating 
the  positions  using  overlapping  plates.  Because  of  the  complexity  here  it  is  helpful 
to  illustrate  the  method  by  a discrete  example.  In  appendix  two  we  give  a numerical 
example  of  two  overlapped  plates  with  ten  stars  using  an  unconstrained  solution.  Here, 
we  show  the  mathematical  constructs  resulting  from  a simple  two-plate  overlap. 

Consider  a “six  constant”  overlap  of  two  plates,  both  having  five  stars  with  three 
of  these  in  common.  There  is  one  reference  star  among  the  three  common  stars;  each 
plate  has  one  other  reference  star.  Altogether,  only  five  stars  will  enter  the  calculations, 
because  the  two  isolated  stars  on  the  plates,  which  are  not  reference  stars,  must  be 
ignored.  We  arrange  the  stars  in  order  of  increasing  right  ascension,  thus  m=5  and 
^=1...5  and  n=2  so  that,  u =1,2.  Figure  2 shows  this  arrangement  of  stars. 

On  plate  1,  the  stars  are  (mi=4);  li=l,  h=2,  ls=3,  l^,=l4=5. 

On  plate  2,  the  stars  are  (a?J2=4);  2i=2,  22=3,  2s  =4,  2m^2m2=24=5 
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There  are  thus  three  reference  stars  ( mr=3);^rj=l,  /<r„^=/^r3=4. 


/ \ 


H = 


\ / 


(^23!  - S$2al 
V2i\  - 
I2jl  - 
V23l  - SV23I 
X2il  — S^2il 
y2tl  - 31)2^1 


— — lilOl 


In  this  case  the  condition  equations  will  be 

J / — S$l,l 

\yi,i  - sj/iii 

XI3I  — S^ljl 

^yiai  - s»/i3i 
*isi  “ *63! 

^yi3i  - ^vi3i 
*1<1  — S^l4l 
Jlil  - SJ/ljl 
2^2,1  — S$2il 

,2/2.1  - SJ/2.1 


— “ijiai 


— — IjlOl 


— -2,1/12 


— —23102 


— -231/12 


— — 24l<l2 


(Oci  - ai)cos^i 

6ci  — ^1 

(ttc3  — as)  cos  63 

^C3  — ^3 

(Oc4  - a4)cOS^4 
^c4  ~ ^4 


First  estimate= 


/ *11 
2/11 
X21 

2/21 

X31 

2/31 

Xil 

yn 

X22 

2/22 

2:32 

2/32 

X42 

2/42 

XS2 

2/52 


\ 


- \ 

- 3T]u 

- 3^21 

- SI/21 

- s6i 

- Sl/31 
~ 0^51 

- SI/M 

- 3^22 

- SI/22 

- 5^32 

- SI/32 

- 0^42 

- SI/42 

- S^52 

- SI/52 
0 

0 
0 
0 
0 
0 


//ill'^ 
(I21 
<^31 
dil 

d22 

d32 

di2 

db2 

0 

0 

0 

0 

0 

Vo/ 


y 


Here,  the  residual  is  d = 


^ ^ / a:  - (a^  + 6r/  + c) 
\y-  {-brj  + ai  + d) 


; optionally  one  would  use 


X — 

y- ST] 

. The  first  matrices  to  examine  are  B,  K,  S and  a 


dimension  of  the  matrix  B is  (16  x 10) 


B = 


/Bai 
0 
0 
0 
0 
0 
0 
0 


0 

B21 

B22 

0 

0 

0 

0 

0 


0 

0 

0 

B31 

B32 

0 

0 

0 


0 

0 

0 

0 

0 

B42 

0 

0 


0 

0 

0 

0 

0 

0 

Bsi 

B52 . 


(3.59) 


. The 


(3.60) 
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We  find  the  individual  elements  in  the  same  way  as  in  the  single  plate  solution.  Using 

W -S' 

V 

equation  (3.24). 


n 241  R — { cos  6 ^ ^ _ 

(xsv,  - d(a,i)  0 ij 


where  W,S,U  and  V are  from 


The  dimension  of  the  matrix  K is  (6  x 10);  K has  the  form 

/I2  0 0 0 o' 


K = 


0 0 I2  0 0 

\ 0 0 0 I2  0 


(3.61) 


The  matrix  S has  dimensions  (16  x 12)  and  is  of  the  form 


■=rT  _ ( “11  “21 

0 


S21  0 

Ssi 

0 

0 

Ssi 

^0  ) 

0 S22 

0 

2s2 

242 

0 

2s2  J 

6)  matrix  = 

Itn' 

1 0 

0 

V 0 

0 

0 ^/IJ/ 

Vttu 

1; 

The  covariance  matrix  is  <7  = 


(To  is  a (6  X 6)  matrix: 


CTx  = 


/<Tll 
0 
0 
0 
0 
0 
0 
\ 0 


0 

(T21 

0 

0 

0 

0 

0 

0 


0 

0 

(722 

0 

0 

0 

0 

0 


0 

0 

0 

(^31 

0 

0 

0 

0 


0 

0 

0 

0 

(732 

0 

0 

0 


O'xx 

0 

0 

0 

0 

0 

0 

(742 

0 

0 


0 

^aa 

0 

0 

0 

0 

0 

0 

(751 

0 


(3.62) 


where  (Tx  is  a (16  x 16)  matrix  and 


0 \ 
0 


0 
0 
0 
0 
0 

(T52  / 


(Ta  = 


(3.63) 


In  this  simple  case  the  matrix  A will  look  as  follows: 


A = 


/B„ 
0 
0 
0 
0 
0 
0 
0 
I2 
0 
\ 0 


0 

B21 

B22 

0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

B31 

B32 

0 

0 

0 

0 

I2 

0 


0 

0 

0 

0 

0 

B42 

0 

0 

0 

0 

I2 


0 

0 

0 

0 

0 

0 

B51 

B52 

0 

0 

0 


“11 

“21 

0 

“31 

0 

0 

S51 

0 

0 

0 

0 


0 \ 
0 

S22 

0 

^32 

S42 

0 

H52 

0 

0 

0 / 


(3.64) 
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First  we  find  the  matrix  L = ^K.  This  is  a (10  x 10)  blcx:k 

diagonal  matrix  which  when  inverted  has  the  form 


L ^ = diag 


s"2(BJi<72i^B21  + BJ2<^22^B22)“^ 
(s^BjjCTgj^  B31  + S“B  J2<732^  B32  + <73)  ^ 
(s^BJ2<742^B42  + <74)"^ 

\s-2(Bjia5-,'B5i  + Bj2aj2'B52)-'  / 


T 


(3.65) 


So  that  BL-^B^  = 


L7i  b7,  0 

0 

0 

0 

0 

0 

0 

0 ^21^22^'21 

B21L22  B^2 

0 

0 

0 

0 

0 

0 H22^22^’21 

B22l'22  B 

0 

0 

0 

0 

0 

0 

0 

0 

B31L33  B3J 

B31L33  B^2  0 

0 

0 

0 

0 

0 

B32L33  B 

B32L33  B J2  0 

0 

0 

0 

0 

0 

0 

0 B42L44BJ2 

0 

0 

0 

0 

0 

0 

0 

0 B51L33  B 

B52 

\ 0 

0 

0 

0 

0 

OB51L-31BJ2 

B52Bj5^B^2  / 
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substituting  into  equation  (3.53)  we  get  J'  = or 


/J'ni 
0 
0 
0 
0 
0 
0 

\ 0 


T' 

•^211 


J', 


’212 


221 

0 

0 

0 

0 

0 


J', 


222 

0 

0 

0 

0 

0 


0 

0 

0 


j; 


0 

0 

0 


’311 


’312 


321 

0 

0 

0 


J', 


322 

0 

0 

0 


0 

0 

0 

0 

0 

I' 

’422 

0 

0 


0 
0 
0 
0 
0 
0 

T'  T' 
•’611  •’612 


Ji 


0 \ 
0 
0 
0 
0 
0 


621 


T' 

•’622 


(3.67) 


where 


'ab 


ab  ’^ab’ 


-InT 


’aa  "ac  ac 


T'  _ 
*^abc 


2.»-1ti  T-lr»T  -1 
ab  ’^abbiaa  ’-’ac^’^ac 


S“£T 


if  b=c 


if  b=c. 


(3.68) 


This  is  a ‘banded  bordered’  matrix  where  the  diagonal  and  adlacent  elemtns  to  the 
diagonal  are  also  filled  on  those  positions  that  correspond  to  stars  imaged  on  more  than 
one  plate. 

The  products  needed  for  the  overlap  solution  are 
S^J'S  = 


(12  X 16)(16  X 16)(16  X 12)  = (12  x 12) 


(3.69) 


(12  X 16)(16  X 16)(16  X 1)  = (12  x 1) 
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Substituting  from  (3.62)  and  (3.67) 

- sV;^bl-^b't<7-^)e 


•z>T  irT 
“11  “21 


0 2?i 


n ■='T 

U _22 


0 2 


WT 

51 


0 ■^T  wT 

U —32  —42 


0 

U “62 


/Jill 

0 

0 

0 

0 

0 

0 

0 \ 

^11 

0 \ 

0 

T' 

•^211 

7' 

J212 

0 

0 

0 

0 

0 

221 

0 

0 

T' 

J221 

7' 

•'222 

0 

0 

0 

0 

0 

0 

222 

0 

0 

0 

7' 

•^311 

7' 

J312 

0 

0 

0 

231 

0 

0 

0 

0 

7' 

J321 

7' 

J322 

0 

0 

0 

0 

232 

0 

0 

0 

0 

0 

7' 

*1422 

0 

0 

0 

242 

0 

0 

0 

0 

0 

0 

7' 

*^511 

Jt>12 

281 

0 

\ 0 

0 

0 

0 

0 

0 

7' 

•'521 

J522  i 

1 0 

252  / 

Hll  IIi2 
Hai  II22 


where 


Hii  - — 7i J'niSii  + Sji J2HS21  + Sji Jaijaai  + 2^j ^ SjjJ^ji2^i 

TT,„ wT  T'  ■z'  I crT  7'  ■=•  I cT  T/  — 

1112  — “21*'211“22  T “31-^312“32  + “61*^612“52  — 

II2I  = 2j2J'221“21  + 232J321S31  + S52J521S51  = ^2  Sj2Jp21* 


fiv 

y~!“ul4ul2^<i2 

HV 


HU 


II22  — “J24222  — 22  + ■='?2J322^32  + J«2^42  + 253 Js22“52  = ^ ^J2^'h22^h2 

fiu 

Also,  STj'd  = ^ 

( “7lJllldll  + “JlJ21ld21  + Sji.T2i2d23  + 2jjJ3jjd3i 

+ “J1J312432  + 2J1  Jjijdsi  + ^Jl^512^S2 
“?2J22ld21  + 2j2J'222d22  + 2^2*^321^31  + 2^24322^32 


+ “42*^422^42  + 252*J52ld51  + 2J2J 


■=*T  T' 


. ■^T  V 


"52*^522'^52 


(3.70) 


(3.71) 


(3.72) 


\ 
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Finally  the  parameters  are  found  from  a = (S^J'S)  where  S^J'S  is  a (12 

X 12)  matrix  and  S^J'd  is  a 12  element  vector.  Therefore  a is  a 12  element  vector, 
i.e.  6 parameters  for  each  plate. 

From  (3.51)  the  equatorial  coordinate  corrections  are  /3  = 5L“^B^cr”^(d  — Ea) 


= S 


y'tii  0 0 0 0 0 0 

0 t21  t22  0 0 0 0 

0 0 0 tsi  t32  0 0 


0 
\ 0 


0 

0 


0 

0 


0 

0 


0 t42  0 

0 0 tgi 


0 \ 
0 
0 
0 

*52  ^ 


( 4ii  — HiiOi  \ 
^21  — H2iai 
^22  — -22«2 
dai  — =3iai 
^32  — “32«2 
^42  — =4202 

4si  — Hsiai 
V4s2  — “5202  / 


where  Uj  = Ljj^Bjaij  therefore: 


13  = s 


( — Elian) 

L22BJiO'2i(421  — “2l02l)  + L22B ^2022(422  ~ “22022) 
L33BT  03i(431  — —3l03l)  + L33BJ2032(432  ~ ^32032) 
L42B  ^2042(442  — E42O42) 

\ L51BJ1  <751(451  — Esiosi)  + L62B ^2^52(452  ~ “52O52)  / 


are  the  ten  weighted  mean  corrections  to  the  star’s  equatorial  coordinates. 


(3.73) 


(3.74) 


Solution  with  Globally  Constrained  Parameters 


This  covers  the  essence  of  the  theoretical  analysis  of  the  overlap  solution.  In 
this  investigation  the  number  of  observations  exceeds  30,000  on  up  to  250  plates  per 
epoch.  This  many  observations  and  plates  will  lead  to  a matrix  that  is  very  difficult  to 
invert,  purely  as  a result  of  the  numerical  computation  loss.  An  alternative  is  to  use 
a different  method  for  finding  the  solution  such  as  Gaussian  elimination.  However,  if 
we  do  not  invert  the  matrix  of  condition  equations  the  covariance  matrix,  and  hence 
the  parameter  errors,  cannot  be  examined.  This  will  restrict  methods  for  checking  the 
solution.  Techniques  need  to  be  examined  to  reduce  the  size  of  the  problem  without 
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an  appreciable  loss  in  accuracy  or  flexibility.  One  method  is  to  restrict  the  number  of 
the  parameters  because  the  matrix  of  condition  equations  in  the  overlap  is  square  with 
dimension  “parameters  times  plates”. 

In  this  investigation  we  exposed  the  McCormick  plates  over  a relatively  short  period, 
one  year  for  the  first  epoch  and  five  months  for  the  second  epoch.  Over  this  short  a 
period  we  can  assume  that  certain  parameters  that  characterize  telescope  properties  will 
remain  constant  from  plate  to  plate.  (This  carries  the  proviso  that  over  the  observation 
period  the  telescope  underwent  no  major  maintenance  or  repairs,  for  example  a cleaning 
of  the  lens.)  This  allows  us  to  restrict  the  overlap  reduction  to  finding  parameters  that 
do  not  vary  from  plate  to  plate  separately  from  those  that  do. 

The  solution  changes  only  in  the  form  of  the  matrix  S.  In  the  general  solution, 
the  V condition  equations  are 


( 

^ \ 

H = 

1 

~ “T'v'nu,*'®" 

/ 

= 0 


where  for  an  eight  constant  model  and  the  i>th  star 


^ 7?  1 0 0 0 ^2  ^7? 

0 0 0 ^ t;  1 ^7/  t; 


n\ 

0 


/“\ 

b 

c 

a' 

b' 

o' 

P 

\q  / 


(3.75) 


(3.76) 


Assume  the  tangential  point  parameters  remain  the  same  across  the  plates,  e.g,  p and  q 
remain  constant  for  all  plates.  In  this  case  the  relationships  (3.75)  will  have  the  form 


H = 


(3.77) 
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where  for  the  uth  star 


^ r?  1 0 0 0 

0 0 0 ^ 7?  1 

/a\ 
b 
c 
a' 
b' 

\c'j 


Ou  = 


-a:  :=) 

-(c:)- 


(3.78) 


Consider  the  form  of  the  matrix  A given  by  (3.36)  the  three  terms  ^ will 

remain  as  given  before.  The  partial  derivative  of  the  condition  equations  H=Q  with 
respect  to  the  plate  parameters,  will  have  a different  form;  in  the  general  case, 


da  5(a,a 


OH, 


/ dH,  dH^  dH^  \ 

bal  dal  ' ' ■ 9<J, 

dH,  dH,  dH, 

da,  da,  ■ ■ ■ da. 


) 


dH.  dHu 

/2.  0 

0 S, 

• • 

• • 

Vo  0 


dHjt, 

dOt/ 

0 \ 

0 


SnJ 


In  the  restricted  case. 


dll 


dH 


da  dia^a^  ■ ■ ■ a„a') 


a//, 

dH, 

dH, 

aoi 

ddg 

da. 

dH^ 

dH^ 

dH, 

da^. 

da^ 

■ ■ ■ da. 

• 

* 

• 

dH„ 

dFf„ 

dH, 

da^ 

da^ 

da. 

0 

...  0 

0 

[I] 

...  0 

* 

1 

I 

\ 0 

0 

■ ■ • “n 

dH,  \ 
da' 
dH, 
da‘ 


dffj^  I 
af  / 


da* 


*2 


(3.79) 


(3.80) 


where  S and  S'  are  as  above. 

The  form  of  the  development  will  be  the  same  as  in  the  generalized  overlap 
discussion.  The  corrections  to  the  reference  stars  will  be  given  by. 


(3.81) 
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and  the  parameters  by 

a = (3.82) 

where  the  matrix  S is  as  in  (3.78). 

This  restricted  approach  can  be  applied  to  many  of  the  parameters.  An  initial  single 
plate  and  overlap  reduction  will  be  carried  out.  Those  parameters  that  have  a standard 
deviation  smaller  than  their  error  will  be  set  constant.  Allowing  these  parameters  to 
vary  adds  nothing  to  the  solution  because  they  vary  less  than  the  noise  of  the  solution, 
and  actually  make  the  solution  worse  because  of  computational  numerical  loss.  A still 
further  extension  of  the  overlap  would  be  to  stochastically  constrain  the  parameters 
that  follow  some  known  distribution  (for  example  a normal  distribution).  This  was 
developed  by  Eichhorn  (1978),  but  it  is  beyond  the  scope  of  this  investigation  to  apply 
these  constraints. 

Example 


Consider  a twelve  constant  overlap  of  the  two  plates,  set  up  as  in  the  unconstrained 
solution.  Make  the  restriction  that  the  tangential  point  corrections  are  the  same  for  both 
plates.  The  derivation  will  proceed  as  before,  with  these  changes 

(3.83) 

where  the  matrix  E is  a (16  x 22)  matrix  of  the  form 


/WT 

“ll 

S21 

0 

[I] 

0 

0 

'Z'T 

“B1 

0 \ 

= 

0 

0 

■Z-T 

“22 

0 

[I] 

[1] 

0 

[1] 

(3.84) 

\ —n 

■Z-T  ' 

wT  ' 

ZtT' 

Z-T  ' 

■Z'T  ' 

[I] 

“21 

“22 

“31 

“32 

“42 

“51 

Each  is  a (2  X 10)  matrix,  thus 


2 -( 
— I 


^nv  1 0 0 0 — 

0 0 0 1 0 


*0  0 (^^1/ + \ 

— nio  IVnu  — T}fiy  / 

(3.85) 


0 
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is  a (2  X 2)  matrix  given  by 

_ f ^nv‘niiv\ 

vlu  ) 

From  (3.51)  the  corrections  to  the  stellar  coordinates  will  be 

/?=  sL-‘B'^<7-^(d-Sa) 


(3.86) 


(3.87) 


where  the  matrix  S is  as  above  and  the  vector  a will  be  of  the  form 

/«\ 

b 

c 
a' 

V 
c' 
e 
/ 

9 

\h/ 


where 


Qu  = 


a = 


For  the  parameters  a — 


(3.88) 


/(lOxl) 

(22  X 1)  = I (10  X 1)  ) and  given  by 
, . . \(2xl) 

(STJ'S)  ^ ((22  X 16)(16  X 16)(16  x 22))"\22  x 16)(16  x 16)(16  x 1)  = 

(22  X 1). 


Using  J from  equation  (3.67),  we  get 


'T 

'11 

S2I 

0 

wT 

>-31 

0 

0 

wT 

“51 

0 \ 

0 

0 

[1] 

0 

zrT 

“32 

[1] 

0 

■Z'T 

“62 

'T  ' 

- 

-Z-T  ' 

wT  ' 

wT  ' , 

'11 

^21 

“22 

“31 

“32 

“42 

“51 

“52  / 

/Jill 

0 

0 

0 

0 

0 

0 

0 

/S„ 

0 

2iA 

0 

T' 

•^211 

7' 

J212 

0 

0 

0 

0 

0 

S21 

0 

[I] 

0 

7' 

•'221 

7' 

J222 

0 

0 

0 

0 

0 

0 

*22 

222^ 

0 

0 

0 

7' 

•^311 

7' 

•'312 

0 

0 

0 

Sai 

0 

^31 

0 

0 

0 

7' 

•'321 

7' 

•'322 

0 

0 

0 

0 

*32 

“32 

0 

0 

0 

0 

0 

7' 

•'422 

0 

0 

0 

*42 

^42"^ 

0 

0 

0 

0 

0 

0 

7' 

•'611 

J612 

^51 

0 

\ 0 

0 

0 

0 

0 

0 

7' 

•'521 

J522  / 

^ 0 

*62 

*52  / 

(3.89) 
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and 


S'^J'S  = 

/ wT  C>T 

“11  “21 

0 n n 

U W37  u u ^51 

0 \ 

0 0 

[1] 

0 wT  -T  Q 

^ “32  “42  ^ 

[>] 

, ■='T'  “T 
\“11  “21 

«— iT  ' * r^T  * <^T  ' * 

•=*22  *='31  ‘='32  ‘='42  '='51 

[1] 

0 

JnjSiT  \ 

J'21i221 

7'  'S' 

J212“22 

7/  “ '_L  7'  'S'  ' 

'^211“21  T *'212“22 

7/ 

•^221^21 

J'522'^22 

J221'^21  + J222'^22 

J311“31 

J312“32 

'Isil'— 31  + Jsi2'^32 

7'  'S' 

J321“31 

Js22'^32 

7/  “ '_L  7'  'S'  ' 

'^321“31  T *'322“32 

0 

*^422'^42 

7/  'S'  ^ 

*^422“42 

'^511'^51 

Jtil2“52 

7/  “ « 1 7/  “ ' 

*'511“61  T>J612“52 

\ J521“81 

7'  'S' 

622“52 

'3c>22'^52  + '^Ml'^Sl  / 

Therefore 


/nn 

II12 

Hi3 

Il21 

II22 

II23 

V II31 

II32 

Has 

where 


n„ 


■Z'T  T/ 


n + 


-r-T  7/ 


“21 


+ 


CT  7/ 
“31*^311 


'31  + 


WT  7/ 


61 


II12 

Has 


rT  7/ 
-21*^212“ 


■=■  j.  wT  1/  -r 
“22  T “3iJ312“32 


wT  -Ti  PE< 
r “61*'512“52 


rT  T/  “ '1  ~T  1/  w ' I “T  II  " ' 1 “T  T/  ™ ' 

■'11"^111“11  + “21*^211“21  T “21'^212“22  + “3lJ311“31 

4-1?T  7/  ■S'„  ' I wT  T/  — > » I “T  7/  “ ' 

' “31*'312“32  T “51'^511“51  + “51‘^612“52 


(3.90) 


(3.91) 


(3.92) 


TT..,  — ’=''^  T'  w -L  -s-T  T'  c?  I ■r-T  7/ 

1121  — “22*'221“21  + “32‘^321“31  + “52’’521 

II22—  “22^222“22  + “32J322“32  + “J2J422“42  + 2^2*^522“52 
H23=  “22(J22iS21  + J222“22)  + 232(  J'32iS3i'+  J322S32') 
"^“42(*^422“420  + “62(3522“62  + J521“5l) 


(3.93) 


51 
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which  can  be  written  as 


3^3'd  = 


^ XT 
V 12 


(3.97) 


With  these  products  we  can  find  the  parameters  and  the  corrections  to  the  stellar 
coordinates. 


Theoretical  Plate  Modelling  Terms 

In  the  single  and  overlap  development  we  use  the  example  of  a six  constant  model. 
In  reality  a more  sophisticated  model  is  needed  to  account  accurately  for  the  complex 
image  formed  by  a telescope.  This  section  discusses  the  various  terms  and  their  physical 
meaning. 

Assume  that  a telescope  has  the  same  projection  properties  as  a pinhole  camera  as 
shown  in  figure  3.  In  this  ideal  system  the  measured  coordinates  (x,y)  with  origin 
at  the  tangential  point  T would  be  related  to  the  standard  coordinates  (^,77)  by  only  a 
factor  of  the  focal  length  of  the  telescope.  In  reality  there  are  many  deviations  from 
the  ideal  situation.  Below  is  a list  of  some  of  the  more  common  ones.  (This  list  is 
not  exhaustive.) 

• axes  rotated  against  the  (x,y)  axes 

• origin  of  the  (x,y)  frame  not  at  the  tangential  point 

• tilt  of  the  ((f,7/)  focal  plane  w.r.t  the  (x,y)  plate  plane  and  noncoincidence  of  the 
optical  axis  with  the  origin  of  the  (^,77)  coordinate  frame 

• magnitude  and  coma  effects 

• radial  distortion 

• tangential  distortion  due  to  components  of  objective  being  improperly  aligned 
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• different  scales  on  each  axis  due  to  different  measuring  screws,  and/or  refraction 
of  the  atmosphere 

• nonorthogonality  of  the  x and  y axes 

• guiding  errors  and  developing  errors  (e.g.,  emulsion  shift) 

• color  terms 

These  deviations  are  not  independent  of  each  other,  and  thus  there  is  a significant 
crossover  of  effects.  The  second  order  effects  are  often  negligible  and  mixed  into 
first  order  effects.  It  is  useful  to  discuss  the  theoretical  form  of  the  more  important 
deviations  we  model  in  this  investigation. 
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Departures  from  a Gnomonic  Projection 
Noncoincidence  of  axes 

The  positioning  of  the  photographic  plate  in  the  telescope  and  then  within  the 
measuring  machine  will  lead  to  a rotation  of  the  plate  with  respect  to  the  frame  of 
reference  on  the  sky.  This  means  the  axis,  defining  the  frame  of  the  sky,  will 
be  rotated  with  respect  to  the  (x,y)  axis,  defining  the  frame  of  the  measuring  machine. 
Noncoincidence  of  axes  A rotation  of  the  ((,r])  frame  will  correct  this: 


— sin  a cos  a 


cos  a sin  q 


(3.98) 


where  s is  the  focal  length. 


Figure  4:  Noncoincidcnce  of  axes 


Origin  shift 


In  the  process  of  transference  from  the  telescope  to  the  measuring  machine  there 
will  be,  in  addition  to  a plate  rotation,  a shift  in  the  origin  of  the  coordinate  frame 
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with  respect  to  the  (x,y)  frame.  If  the  origin  of  the  (^,7/)  frame  in  the  (x,y)  frame  is 

(c,cO,  then  this  can  be  allowed  by  a constant  subtraction  from  the  (x,y)  coordinates. 

Combining  this  change  into  equation  (3.98)  yields; 

f X — c\  ( cos  a sin  Q \ / f \ 

, =4  . r . (3.99) 

\ 2/  ~ \ — sin  a cos  Q J \ r}  J 


11 

I 

y 

A 


(c,  c' ) 


X-* 


Figure  5:  Origin  Shift 


Incorrect  scale 


In  an  ideal  system  the  (^,7/)  coordinates  would  be  related  to  the  (x,y)  coordinates  by 
the  focal  length;  i.e.  (x,y)  = s (^,7;).  In  practice  the  actual  focal  length  of  the  telescope 
will  vary  due  to  temperature  changes  and  everyday  use.  We  therefore  treat  this  as  a 
variable  and  correct  by  incorporating  a small  modification  in  the  direct  ^ and  7;  terms. 
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The  relationship  between  standard  and  rectangular  coordinates  now  has  the  form 


/ X — c / cos  a sina\/^\  f ^ \ f 

\y  — c'J  V — sina  cos  a J\i]J  \a'J\T)J 


(3.100) 


where  a and  a'  are  constants. 

These  terms  can  be  combined  if  we  include  the  cosine  and  sine  terms  as  constants, 
then  the  deviations;  noncoincidence,  origin  shift  and  scale  difference  in  matrix  format 


are 


/ X — c\  / cos  a sina\  / f ^ \ f 

\y  ~ J V~sina  cos  a J\rjJ  \a' J \i]  J 


= s 


cos  a + SI]  sin  a — + c 

sin  a + ST]  cos  a — a'r]  + c' 
(coso  - + (sina);?  + f 

cos  a — y j T]  — (sin  a)^  + y 
/a\ 


= s 


b 

c 

a' 

b' 

\c'J 


(3.101) 


where  a = cos  a — 6 = sin  q c = -,  a'  = — sin  a,  b'  = cos  a — — , and  c'  = — . 

s’  s’  ’ s’  s 

This  is  the  six  constant  model.  The  rest  of  this  discussion  will  build  on  this  model; 


it  is  therefore  important  to  first  say  a word  about  this  idea  of  modelling.  From  equation 
(3.101)  it  would  appear  that  if  our  projection  were  only  subject  to  the  deviations:  non- 
conjunction of  the  axis  and  scale  differences,  then  b=-a'.  However  this  does  not  happen 
even  in  this  simple  case.  This  is  because  the  above  analysis  hides  second  order  effects, 
and  we  have  made  no  provisions  to  account  for  atmospheric  refraction.  There  will  be 
second  order  mixing  and  the  reader  should  be  aware  that  these  physical  interpretations 
of  the  models  are  only  good  to  the  first  order.  This  will  be  alleviated  as  we  add  more 
terms  that  are  more  difficult  to  model. 
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Focal  plane  tilt  and  an  incorrect  tangential  point 


If  either  the  photographic  plate  is  tilted  with  respect  to  the  focal  plane  or  the  assumed 

tangential  point  is  not  the  same  as  the  intersection  of  the  optical  axis  with  the  plate,  the 

differential  change  in  (x,y)  with  respect  to  the  plate  center  (oo.^o)  is 

d{x,y)  d{x,y)  d{^,n) 


d(aoJo)  d{^,T})  d{ao,So) 
d(x,y)  _ ^ a b 
b', 

d(Qo,6o)  5(E,H,Z)  5(ao,^o)  ' 


From  equation  (3.2)  we  get 


^ Z'  ^ Z 


dU,y)  ^1/1  0 -n 
(H,H,Z)  ZVO  1 -vJ 


and 


T = 


/H\ 


H 


\zy 

5(ao,6o) 


9(H,H,Z) 
9(ao,  io) 


This  simplifies  to 


d{ 


cos  6 sin  (a  — Qo) 

cos  5o  sin  6 — sin  So  cos  S cos  (a  — oo ) 

V sin  6o  sin  6 + cos  So  cos  S cos  (a  — ao)  / 

— cos  S cos  (a  — ao ) 0 

— sin  So  cos  S sin  (a  — ao)  — sin  So  sin  S — cos  6o  cos  S cos  (a  — ao) 
cos^o  cos  # sin  (a  — ao)  sin  So  cos  S — sin^o  cos  6 cos  (a  — ao) 

' H sin  So  — Z cos  So  0 \ 

— E sin  So  Z . 

H cos  So  H / 


di^,v)  ^ d{^,r))  d(E,li,Z) 

d{ao,6o)  5(E,n,Z)  c»(ao,<5o) 

/ H sin  So  — Z cos  So  0 \ 
— Esin^o  Z 


j_n  0 
1 


V 


E cos  So 

T]  sin  ^0  - cos  So  - ^ ^ cos  So 
— ^sin^o  — ^1]  cos  So 


H/ 


-ri.) 


(3.102) 


(3.103) 


(3.104) 


(3.105) 
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1 and  2 have  positions  that  depend 
both  on  their  ^ and  r|  values. 


Figure  6:  Effect  of  Plate  Tilt  or  Incorrect  Tangential  Point 


and  finally, 

d(x,y)  6 \ /Tjsinio  - cos6o -^^cos(5o  \ 

^(ooi^o)  Va'  6 / \ -^sm6o  — ^ycosSo  — 

( a{r)s\i\6o  — cos^o  - cos6q)  — b{—^  sin^o  - Ol^os6o)\ 

/ —a  cos  So  + (6  sin  5q)^  + (a  sin  ^o)*?  + (bcos6o)^ri  — (a  cosfioX^  '\ 

"V  b'-a'^  + b'rf  ) 

(cf.  Eichhorn,  1974). 

Most  models  will  at  a minimum  contain  the  linear  six  constant  model.  (If  we 
have  an  ideal  measuring  machine  then  we  can  use  a four  constant  model  if  the  effects 
of  differential  atmospheric  refraction  are  removed  beforehand.)  The  constant  term 
— acos^o  and  the  t]  term  asin^o  in  (3.106)  can  be  incorporated  in  the  six  constant 
c'  and  b parameters  respectively.  The  only  new  terms  that  need  to  be  introduced  are 
those  that  are  functions  of  and  ?/“. 

These  terms,  introduced  by  a tangential  point  shift,  are  also  those  terms  that  occur 
in  plate  (focal  plane)  tilt.  We  illustrate  this  by  looking  at  the  geometry.  Tilting  of  the 
plate  causes  the  coordinate  system  to  have  different  scales  depending  on  the  position 
in  the  plate. 
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The  ^x.?/x  coordinates  of  the  point  x will  be  functions  of  the  ideal  (gnomonic)  ( 
and  rj  respectively,  e.g  ^x=<fx(6^)  /;x=?/x(6^/)-  So  the  tilt  will  add  square  and  cross 


where  p and  q are  constants. 
Magnitude  and  coma  effects 


The  stellar  magnitude  terms  and  lens  coma  will  cause  the  image  of  a star  to  shift. 
This  is  particularly  important  in  studies  that  use  concentric  plates  as  these  terms  will 
cause  an  expansion  of  the  observed  region.  The  reason  is  simple,  a brighter  star  will 
cause  the  centroid  found  for  an  image  to  move  out  along  a radial  vector.  Figure  7 
shows,  in  a highly  exaggerated  example,  this  effect  for  two  stars;  the  radial  vector  is 
pointing  down. 

The  two  images  should  have  the  same  position  for  their  centroid,  but  because  of 
coma  the  best  fit  ellipse  in  the  brighter  star  is  centered  further  in  along  the  radial  vector 
than  that  for  the  dimmer  star.  We  model  these  effects  by  including  the  terms 


terms  in  ^ and  r;  into  the  model.  So  we  can  allow  for  both  these  effects  by  introducing 


the  terms 


(3.107) 


111  — (m) 
0 


0 (m-{m))i  (m-(m})^2 

m-{m)  {m  — {m))rj  {m-{m))rf 


/e\ 

/ 


(3.108) 


9 


\h) 


where  e,  f,  g,  and  h are  constants. 
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Figure  7:  Effect  of  Magnitude  and  Coma  Terms 


Radial  distortion 


Barrel  and  pincushion  distortion  are  lens  aberrations  which  produce  radial  terms. 
To  correct  for  this  distortion  a term  of  the  form 


Ax 


Ay 


\riii^  + t) 


(0 


(3.109) 


where  / is  a constant,  is  included  in  many  models. 

This  completes  the  discussion  of  the  more  important  departures  from  a gnomonic 
projection  that  will  affect  the  position  of  a star.  There  are  many  other  possible  systematic 
departures  from  a gnomonic  projection,  color  terms,  magnitude  squared  terms,  radial 
magnitude  terms  to  list  a few.  Analysis  of  the  residuals  will  highlight  those  terms  to 


include. 
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Typical  Models 


Assume  the  relationship  of  measured  coordinates  to  standard  coordinates  is  of  the 


form 


y 


= s\  ^ j + Sa,  where  S is  the  model  matrix  and  a is  the  parameter  vector. 


Below  we  examine  the  6,  8,  12,  and  16  constant  models.  Other  models  are  also  tested 
in  this  study,  but  the  ones  listed  are  good  examples  of  the  combinations  tried. 


The  6 constant  model 


We  have  already  discussed  the  6 constant  model  and  its  physical  interpretation. 


r;  1 0 0 ON 

Vo  0 0 ^ ij’ 


6 


(3.110) 


• c and  c'  account  for  differences  between  the  origins  of  the  measurable  rectangular 
image  coordinates  (x,y)  and  the  standard  coordinates  (^,r;). 

• a and  b'  allow  for  an  incorrect  focal  length  and  different  x and  y scales  as  caused 
by  differential  refraction. 

• b and  a'  allow  for  rotation  of  the  x,y  frame  vs  the  frame  and  non-perpendicularity 


of  the  axes. 
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The  8 constant  model 


For  the  8 constant  model  we  just  add  linear  terms  in  the  magnitude. 


/^»jl000  »n  — (m)  0 \ 

\000^»7l  0 m—  (m)  ) ’ 


(a\ 

b 


c 

a! 


V// 


(3.111) 


• e and  f allow  for  magnitude  effects 

This  model  is  only  really  important  to  use  as  a base  for  introducing  global  param- 
eters. It  corrects  for  the  major  plate-to-plate  dependent  terms  resulting  from  guiding, 
scale  and  plate  placement  errors.  Additional  parameters  may  be  considered  global,  and 
should  be  tested  for  each  set  of  plates  for  that  possibility. 

The  12  constant  model 


If  we  add  in  terms  to  correct  for  plate  tilt,  coma  and,  radial  effects  we  get  the 
12  constant  model  as  shown  in  (3.112).  In  this  model  the  tilt  and  tangential  point 
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corrections  are  assumed  to  have  dependencies  that  mirror  in  the  x and  y directions. 


n r,  i 0 0 0 e 

\0  0 0^  f]  I ^T)  i]~ 


m — (m) 
0 


0 

m — (m) 


/ a \ 
6 
c 
a' 
b' 
c' 

P 

9 

e 

/ 

9 

\hj 


(m  - {m))rj  + Tf^)  ) 


(3.112) 


• p and  q allow  for  tilt  and  tangential  point  corrections 

• e and  f allow  for  magnitude  effects 

• g corrects  for  coma,  it  should  have  the  same  coefficient  in  both  x an  y 

• h corrects  for  radial  distonion,  again  the  coefficient  should  be  the  same. 

The  16  constant  model 


The  most  sophisticated  model  usually  attempted  is  the  sixteen  constant  model.  Here 
we  allow  the  tilt  and  tangential  point  corrections  to  have  different  x and  y dependencies 
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(parameters). 
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(3.113) 


* r,  p',  q',  and  r'  allow  for  separate  tilt  and  tangential  point  corrections. 


This  model  proved  to  be  dominated  by  the  cubic  radial  term  and  another  form  of  it 

replaces  this  term  with  a second  order  magnitude  modelling  term,  ( ~ ^(h) 

\{m  — {m))T)‘‘  J ' ’ 


Summary  of  potential  model  terms 


The  terms  that  might  be  included  can  be  simply  split  into  three  categories.  Terms 
that  usually  require  different  parameters  in  the  functional  forms  of  the  x and  y coor- 
dinates; terms  that  may  have  the  same  parameter  in  both  the  x and  y values,  but  may 
also  be  split;  and  finally  terms  that  usually  have  the  same  parameters  in  x and  y,  i.e. 
are  radial  in  nature,  and  will  not  usually  be  split. 
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Terms  that  usually  have  different  x and  y parameters  include  all  the  linear  terms 
in  the  standard  coordinates  and  the  magnitude:  ?/,  and,  (m-(m)).  These  are 

strongly  dependent  upon  the  conditions  of  that  particular  exposure,  e.g.  centering  of 
the  plate,  guiding  during  exposure,  and  scale  variations  (which  are  in  turn  dependent 
on  the  temperature  at  the  time  of  exposure).  These  terms  are  then  usually  found  for 
each  plate  and  each  coordinate  on  those  plates;  it  would  be  unwise  to  try  and  globally 
fit  them  or  constrain  them  to  the  same  value  for  both  coordinates.  Because  of  their 
usually  small  random  variation  the  scale  terms,  the  direct  terms  in  ( and  rj,  are  very 
good  candidates  for  stochastic  constraints  because  they  would  usually  be  assumed  to 
be  normally  distributed  around  some  mean  value. 

Model  terms  that  are  sometimes  different  in  the  x and  y parameters  and  sometimes 
constrained  to  the  same  value  for  both  the  coordinates  are  the  second  order  terms 
7/^,  and  ^77)  and  — if  they  are  included  — cubic  terms  ((^^,  7/^  and,  These 
allow  for  plate  tilt,  tangential  point  corrections,  some  distortion  effects  and  differential 
atmospheric  refraction;  they  are  usually  small.  If  an  analysis  of  their  values  in  a single 
plate  solution  shows  that  they  vary  widely  for  the  two  coordinates,  then  they  can  be  split. 
These  parameters  are  very  telescope  dependent,  and  providing  the  telescope  undergoes 
no  overhauls  during  the  period  of  exposure  for  a set  of  plates,  these  parameters  will  be 
good  candidates  for  global  constraints. 

The  distortion  terms  and  and  coma  terms  [(m  — (ttz))^, 

{m  — {m))rj,  {m  — {m  — {m))r]“,  and  {m  — (777))^]  can  usually  be  assumed 

to  be  the  same  in  both  the  x and  y coordinate.  These  terms  are  mainly  related  to  the 
lens,  and  they  may  be  considered  either  globally  constrained,  stochastically  constrained 
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or  removed  in  some  pre-overlap  process.  The  nature  of  the  actual  problem  indicates 
the  road  to  follow. 

As  will  be  discussed  in  the  next  chapter  the  use  of  many  parameters  may  introduce 
large  variances  in  the  possible  solution,  while  leaving  out  terms  introduces  the  risk  of 
not  including  effects  that  may  be  very  important  to  model.  An  analysis  of  parameter 
variance  by  Eichhorn  and  Williams  (1963)  shows  that  including  unneeded  terms  will 
increase  the  ultimate  error  of  the  stellar  position  by  the  introduction  of  errors  through 
the  parameter  variances. 


CHAPTER  4 
OBSERVATIONS 


We  have  obtained  observations  from  three  epochs.  The  first  epoch  (the  Astrographic 
Catalogue,  AC)  is  available  in  the  form  of  published  x and  y coordinates.  The  second 
and  third  epochs  required  measuring  of  photographic  plates  on  an  automatic  plate 
measuring  machine.  The  first  two  sections  in  this  chapter  will  discuss  the  observations, 
and  the  last  section  will  discuss  the  measuring  of  the  plates. 

Collection  of  the  Observations 

The  First  Epoch 

For  the  first  epoch  we  used  the  Astrographic  Catalogue  (AC).  The  AC  was  one  of 
the  largest  astronomical  enterprises  ever  undertaken.  In  1890  twenty-two  observatories 
all  over  the  world  started  a program  to  map  the  sky  using  a ‘standard’  astrograph.  A 
standard  astrograph  is  a double  telescope  whose  photographically  corrected  lens  has  a 
plate  scale  of  r=lmm.  The  hope  was  that  by  requiring  this  adherence  to  a particular 
instrument  and  by  imposing  strict  guidelines  on  the  measurements  of  the  stellar  image 
coordinates  there  would  be  a consistent  map  of  the  whole  sky  for  an  epoch  around  1900. 

The  AC  has  provided  an  invaluable  reference  source  for  astronomers  ever  since.  In 
particular  it  provides  a long  baseline  and  has  a strong  corner-in-center  overlapping  plate 
pattern.  The  observation  of  the  sky  was  finally  completed  in  1950,  and  the  measuring 
soon  after. 

The  fruits  of  the  whole  enterprise  were  published  in  the  form  of  a multi-volume 
catalogue  with  plate  parameters  for  a six  constant  model  and  the  x and  y coordinate 
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and  magnitude  for  each  star  on  that  particular  plate.  The  user  could  then  from,  the 
original  measurements  to  determine  the  stars’  equatorial  coordinates  and  possibly  further 
improve  upon  if  then  needed.  The  catalogue  has  a limiting  magnitude  between  12-1 3th 
magnitude,  although  with  today’s  measuring  devices  this  magnitude  limit  could  probably 
be  extended. 

The  United  States  Naval  Observatory  is  currently  entering  all  the  coordinates  and 
attempting  to  re-reduce  the  entire  catalogue  with  a block  adjustment.  The  author  is 
in  debt  to  the  US  Naval  Observatory  for  providing  the  AC  of  the  Orion  region  from 
their  database.  The  reference  catalogue  used  is  the  Astrographic  Catalogue  of  Reference 
Stars  (ACRS),  also  provided  by  the  US  Naval  Observatory;  who  produced  this  catalogue 
solely  for  the  purpose  of  re-reducing  the  whole  of  the  Astrographic  Catalogue.  The 
density  of  the  ACRS  is  an  average  of  8 stars  per  square  degree,  so  nominally  there 
will  be  32  on  each  Astrographic  plate  and  8 on  each  McCormick  plate  (because  of  the 
bright  star  density  in  Orion  this  density  is  actually  higher  for  this  study). 

Figure  8 shows  the  pattern  of  the  San  Fernando  and  Algiers  photographic  plates. 
The  three  belt  stars  of  Orion  are  shown  for  orientation.  Note  the  strong  overlapping 
comer  in  center  pattern.  This  means  that  a star’s  image  can  be  on  up  to  5 separate 
plates. 

The  San  Fernando  region 

We  used  19  plates  from  the  San  Fernando  region.  Table  1 summarizes  the  original 
AC  data.  Listed  are  a running  plate  number,  astrographic  catalogue  plate  number,  right 
ascension  in  hours,  declination  in  degrees,  epoch  of  Observation  and  the  six  parameters 
given  in  the  AC. 
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Figure  8:  First  Epoch  Plate  Orientation 


By  examining  the  epochs  of  the  plates  one  can  see  that  the  majority  of  them  were 


taken  around  1893,  and  then  there  is  another  cluster  of  observations  in  1909.  The 


latter  plates  only  cover  stars  above  0°  declination;  this  will  be  drawn  upon  in  the  final 


summary.  For  now,  we  shall  just  note  that  the  San  Fernando  region  was  made  up  of 
material  clustered  around  two  epochs,  1893  and  1909. 


The  first  and  fifth  parameters,  a and  e,  are  quite  large  and  positive;  this  implies  that 
the  focal  length  of  the  San  Fernando  astrograph  was  approximately  37mm  smaller  than 
the  standard  astrographic  focal  length  of  3437mm. 
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Table  1:  San-Femando  Aslrographic  Plates  Data 


# AC#  Qph  Sq°  Epoch  Six  constant  model  published  with  the  AC 


1 

2286 

5.33 

2.0 

1910.0 

37.75 

3.02 

0.144 

-3.20 

37.29 

0.105 

2 

2287 

5.47 

2.0 

1910.0 

37.62 

9.69 

0.028 

-9.94 

37.17 

0.068 

3 

1954 

5.27 

1.0 

1909.0 

37.22 

-12 . 5 

-.552 

12.33 

36.49 

-.036 

4 

1956 

5.40 

1.0 

1909.0 

37.64 

0.80 

-.445 

-0.87 

36.88 

-.124 

5 

1957 

5.53 

1.0 

1909.0 

37.20 

9.13 

0.022 

-9.24 

36.44 

-.018 

6 

0621 

5.20 

0.0 

1894.0 

37.45 

-8.72 

0.454 

8.75 

36.89 

-1.35 

7 

0135 

5.33 

0.0 

1892.1 

38.79 

-1.70 

-.948 

1.70 

38.27 

0.179 

8 

0630 

5.47 

0.0 

1894.0 

36.80 

6.77 

0.242 

-6.77 

36.28 

-.308 

9 

0136 

5.60 

0.0 

1892.1 

37.66 

-13.6 

-.319 

13.62 

37.11 

-.409 

10 

0411 

5.27 

-1.0 

1893.1 

37.36 

-42.6 

1.023 

42.80 

36.77 

-.515 

11 

0412 

5.40 

-1.0 

1893.1 

37.33 

4.38 

1.455 

-4.31 

36.77 

-.475 

12 

0080 

5.53 

-1.0 

1892.0 

37.19 

-19.5 

-.139 

19.56 

36.60 

0.099 

13 

0413 

5.67 

-1.0 

1893.1 

37.27 

-10.4 

0.202 

10.52 

36.68 

-.263 

14 

0403 

5.80 

-1.0 

1893.1 

36.58 

-40.7 

0.272 

40.59 

35.99 

-.589 

15 

0424 

5.20 

-2.0 

1893.1 

37.13 

-18.2 

0.418 

18.10 

36.54 

-.483 

16 

0400 

5.33 

-2.0 

1893.0 

36.69 

7.36 

-.297 

-7.36 

36.07 

1.530 

17 

1108 

5.47 

-2.0 

1896.0 

37.43 

2.36 

0.093 

-2.33 

36.84 

-.897 

18 

0161 

5.60 

-2.0 

1892.1 

37.84 

-7.16 

-.289 

7.19 

37.21 

-.442 

19 

0163 

5.73 

-2.0 

1892.1 

37.01 

-6.32 

0.297 

6.18 

36.38 

-.196 

The  second  and  forth  parameters,  b and  d,  vary  quite  widely  and  are  usually  the 
opposite  of  each  other.  These  represent  the  rotation  terms  and  this  variation  is  expected 
if  one  also  takes  into  account  that  the  third  and  sixth  parameters  remain  quite  small. 
The  San  Fernando  observers  and  measurers  were  very  careful  to  ensure  the  plate  was 
well  centered  and  reasonably  careful  to  ensure  consistent  alignment  of  the  plates  with 
respect  to  the  sky. 
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The  Algiers  region 


From  the  Algiers  observatory  we  used  12  plates.  Table  2 summarizes  the  original 
AC  data  listed  in  the  same  format  as  table  1. 

Table  2:  Algiers  Astrographic  Plates  Data 
# AC#  Ooh  6o°  Epoch  Six  constant  model  published  with  the  AC 


20 

2450 

5.27 

-3.0 

1896.1 

-1.00 

2.68 

0.608 

-2.75 

-1.65 

0.128 

21 

4223 

5.40 

-3.0 

1909.1 

-1.10 

0.69 

0.062 

-1.13 

-1.68 

-.212 

22 

4224 

5.53 

-3.0 

1909.1 

-1.03 

0.89 

0.881 

-1.31 

-1.89 

-.107 

23 

2453 

5.67 

-3.0 

1896.2 

-1.13 

0.10 

0.033 

-0.38 

-1.96 

0.131 

24 

1498 

5.80 

-3.0 

1894.0 

-1.51 

-1.24 

0.183 

1.34 

-2.06 

-.145 

25 

1697 

5.33 

-4.0 

1894.2 

-1.24 

-1.65 

0.141 

1.48 

-1.55 

0.004 

26 

2451 

5.47 

-4.0 

1896.2 

-0.41 

-1.06 

-.535 

0.72 

-1.86 

2.729 

27 

1716 

5.60 

-4.0 

1894.2 

-1.00 

-0.96 

0.079 

0.10 

-1.82 

-.036 

28 

1728 

5.73 

-4.0 

1894.2 

-0.83 

-2.71 

-.138 

2.30 

-0.96 

-.091 

29 

2429 

5.40 

-5.0 

1896.1 

-0.62 

-0.31 

2.604 

-0.07 

-2.13 

-1.38 

30 

1496 

5.53 

-5.0 

1894.0 

-1.48 

1.00 

0.216 

-1.03 

-1.99 

0.298 

31 

1512 

5.67 

-5.0 

1894.0 

-1.17 

-2.20 

0.268 

1.82 

-1.96 

-.167 

Once  again  the  observations  are  clustered  around  two  epochs,  1895  and  1909. 
For  this  study  we  must  pick  one  epoch  as  a mean  epoch  of  observation  and  ignore 
the  fact  that  the  plates  were  actually  taken  at  separate  epochs.  This  is  not  strictly 
necessary;  because  we  have  three  epochs  of  observation,  it  would  be  possible  to  include 
a proper  motion  term  with  each  observation.  This  would  however  vastly  increase  the 
computational  work  required  and  the  number  of  free  parameters  per  image.  Considering 
that  the  literature  cites  the  nominal  measuring  error  of  the  AC  at  four  microns,  the 
minimum  proper  motion  to  be  measurable  above  the  measuring  error  would  be  0.02 
arcseconds  per  year.  There  are  very  few  stars  with  this  high  a proper  motion,  so  we 
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will  use  the  simple  approach  of  assuming  the  observations  to  be  made  at  the  middle 
epoch,  which  for  these  data  is  approximately  1900. 

The  first  and  fifth  parameters,  a and  e,  are  small;  this  implies  that  the  focal  length  of 
the  Algiers  astrograph  was  close  to  the  standard  astrographic  focal  length  of  3437mm. 

The  second  and  forth  parameters,  b and  d,  do  not  vary  much  and  also  the  third  and 
sixth  parameters  are  small.  The  Algiers  observers  and  measurers  were  very  careful  in 
both  centering  and  orienting  the  photographic  plates  during  observing  and  measuring. 

The  Second  Epoch 

For  the  second  epoch  we  used  photographic  plates  exposed  by  Heinrich  Eichhom  at 
the  26-inch  refractor  of  the  Leander  Me  Cormick  Observatory.  This  is  an  exceptional 
refractor  with  a long  history  in  parallax  determinations  and  other  astrometric  research. 
With  a plate  scale  of  227mm  it  has  a much  larger  plate  scale  than  the  standard  astrograph 
giving  more  precise  information  per  unit  area. 

During  these  observations  Eichhorn  exposed  each  region  four  times,  rotating  the 
plate  once  after  the  second  exposure.  This  was  for  a number  of  reasons:  more  images  of 
each  star  would  reduce  the  centering  error  that  is  inherent  in  any  measuring  process;  by 
rotating  the  plate  he  nullified  the  effects  of  emulsion  slip,  and  any  systematic  observing 
errors  that  depended  on  the  orientation;  and  finally  he  ensured  that  if  any  particular 
exposure  was  subject  to  a large  guiding  error  there  would  be  other  exposures  of  that 
region.  Figure  9 illustrates  the  orientation  of  the  67  regions  he  choose  to  observe. 

Again  note  the  overlapping  pattern:  each  star  will  appear  in  a maximum  of  four 
regions,  and  because  of  the  repeated  exposures,  each  star  has  the  potential  to  be 
measured  up  to  16  times.  This  gives  a strong  tie  for  the  overlapping  method  to  exploit. 
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Figure  9:  Second  Epoch  Plate  Orientation 


Table  3 summarizes  the  region  data.  Listed  are  running  plate  number,  McCormick 
Observatory  plate  number;  right  ascension  in  hours;  declination  in  degrees;  date;  eastern 
standard  time;  hour  angle;  seeing  (l=high  seeing,  5=low  seeing);  transmittance  or  cloud 
cover  (l=very  covered,  5=clear);  and  the  temperature  of  the  lens  at  the  end  of  the 


exposure. 


MC# 

66807 

66810 

66812 

66814 

66816 

66817 

66833 

66835 

66837 

66839 

66841 

66843 

66853 

66855 

66857 

66863 

66865 

66884 

66886 

66895 

66897 

66899 

66901 

66911 

66913 

66915 

66917 

66919 

66921 

66923 

66925 

66927 

66929 

66932 

66933 

66942 

66944 

66946 

66948 

66957 
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Table  3:  Second  Epoch  McCormick  Observatory  Plates  Data 


Ooh 

So^ 

Epoch 

EST 

HA 

s 

T 

Temperature 

5.66 

-3.6 

1955.7 

09:13 

0 : 40E 

4 

4 

36 

5.61 

-3.7 

1955.7 

09:39 

0:11E 

4 

4 

36 

5.55 

-3.6 

1955.7 

10:01 

0:17W 

4 

4 

36 

5.55 

-3.1 

1955.7 

10:22 

0:36W 

4 

4 

36 

5.61 

-3.1 

1955.7 

10:45 

0:54W 

4 

4 

36 

5.65 

-3.1 

1955.7 

11:00 

1:06W 

4 

4 

36 

5.40 

-3.6 

1956.1 

08:28 

0:56E 

5 

4 

35 

5.46 

-3.6 

1956.1 

08:48 

0:39E 

5 

4 

35 

5.50 

-3.7 

1956.1 

09:19 

0:10E 

5 

3 

35 

5.50 

-3.1 

1956.1 

09:51 

0:22W 

5 

3 

35 

5.46 

-3.1 

1956.1 

10:23 

0:56W 

5 

4 

35 

5.71 

-3.1 

1956.1 

10:51 

1 : 09W 

5 

4 

35 

5.40 

-3.1 

1956.1 

08:27 

0:51E 

5 

4 

36 

5.36 

-3.1 

1956.1 

08:47 

0:28E 

5 

4 

36 

5.35 

-2.6 

1956.1 

09:10 

0:05E 

5 

4 

36 

5.51 

-2.6 

1956.1 

10:19 

0 : 54W 

5 

4 

36 

5.71 

-2.6 

1956.1 

10:48 

1 : 12W 

5 

4 

36 

5.26 

-2.1 

1956.1 

07:53 

0:56E 

4 

4 

38 

5.26 

-1.6 

1956.1 

08:01 

0 : 48E 

4 

4 

38 

5.41 

-2.6 

1956.1 

08:16 

0:15E 

4 

5 

44 

5.46 

-2.6 

1956.1 

08:34 

0:00E 

4 

5 

44 

5.56 

-2.7 

1956.1 

08:55 

0:16W 

4 

5 

44 

5.61 

-2.6 

1956.1 

09:11 

0:28W 

4 

5 

44 

5.31 

-2.6 

1956.1 

07:40 

0:40E 

5 

4 

48 

5.31 

-2.1 

1956.1 

07:58 

0:22E 

5 

4 

48 

5.36 

-2.1 

1956.1 

08:29 

0 : 05W 

4 

4 

48 

5.40 

-2.1 

1956.1 

08:45 

0 : 18W 

4 

3 

48 

5.30 

-1.6 

1956.1 

07:05 

1:11E 

4 

4 

50 

5.36 

-1.6 

1956.1 

07:25 

0:54E 

4 

4 

50 

5.40 

-1.6 

1956.1 

07:46 

0:36E 

4 

4 

50 

5.45 

-1.6 

1956.1 

08:07 

0:18E 

4 

4 

50 

5.50 

-1.6 

1956.1 

08:29 

0 : OIW 

3 

4 

50 

5.55 

-1.6 

1956.1 

08:48 

0:17W 

3 

4 

50 

5.55 

-2.1 

1956.1 

09:12 

0 : 41W 

2 

3 

50 

5.61 

-2.2 

1956.1 

09:34 

0 : 59W 

2 

3 

50 

5.24 

-1.1 

1956.1 

06:56 

1 : 06E 

5 

5 

45 

5.30 

-1.1 

1956.1 

07:13 

0:51E 

5 

5 

45 

5.35 

-1.1 

1956.1 

07:35 

0:32E 

5 

5 

45 

5.41 

-1.1 

1956.1 

07:53 

0: 17E 

5 

5 

45 

5.45 

-2.1 

1956.1 

07:26 

0:43E 

4 

5 

42 

75 


# 

MC# 

Ooh 

o 

o 

Epoch 

EST 

HA 

s 

T 

Temperature 

41 

66959 

5.50 

-2.1 

1956.1 

07:44 

0:28E 

4 

5 

42 

42 

66963 

5.35 

-0.1 

1956.1 

06:51 

0:52E 

3 

5 

52 

43 

66965 

5.41 

-0.1 

1956.1 

07:08 

0:38E 

5 

5 

52 

44 

66967 

5.45 

-0.1 

1956.1 

07:28 

0 : 21E 

5 

5 

52 

45 

66969 

5.45 

-0.6 

1956.1 

07:48 

0:01E 

5 

5 

52 

46 

66971 

5.51 

-1.1 

1956.1 

08:06 

0:17W 

5 

5 

52 

47 

66976 

5.56 

-1.1 

1956.1 

08:51 

0:56W 

5 

5 

52 

48 

66978 

5.61 

-1.6 

1956.1 

09:15 

1:16W 

5 

5 

52 

49 

66979 

5.36 

0.5 

1956.1 

06:48 

0:47E 

3 

4 

44 

50 

66981 

5.41 

0.5 

1956.1 

07:05 

0:33E 

3 

4 

44 

51 

66983 

5.45 

0.5 

1956.1 

07:25 

0:16E 

3 

4 

44 

52 

66985 

5.51 

0.4 

1956.1 

07:43 

0:01E 

3 

4 

44 

53 

66988 

5.55 

-0.1 

1956.1 

08:14 

0:27W 

3 

3 

44 

54 

66990 

5.60 

-1.2 

1956.1 

08:36 

0 : 45W 

3 

3 

44 

55 

66991 

5.51 

-0.1 

1956.1 

06:47 

0-.46E 

4 

4 

42 

56 

66994 

5.51 

-0.6 

1956.1 

07:04 

0:29E 

4 

4 

39 

57 

66996 

5.55 

-0.6 

1956.1 

07:29 

0:07E 

4 

4 

39 

58 

66997 

5.61 

-0.6 

1956.1 

07:41 

0:  OIW 

4 

4 

39 

59 

66999 

5.65 

-1 . 6 

1956.1 

08:05 

0:08W 

4 

4 

39 

60 

67001 

5.70 

-2.2 

1956.1 

08:32 

0:46W 

5 

4 

37 

61 

67003 

5.70 

-1.7 

1956.1 

62 

67005 

5.35 

-0.6 

1956.1 

63 

67007 

5.41 

-0.6 

1956.1 

64 

67008 

5.35 

0.9 

1956.1 

65 

67010 

5.41 

0.9 

1956.1 

66 

67012 

5.46 

0.9 

1956.1 

67 

67015 

5.51 

0.9 

1956.1 

Figure  10  is  a picture  of  the  refractor.  The  photographic  plate  holder  is  shown  quite 
clearly  at  the  base  along  with  the  offset  guiding  mechanism  and  exposure  timer.  At  the 
top  of  the  telescope  one  can  see  the  housing  for  both  the  lens  cover  and  the  objective 
grating  used  in  the  second  and  third  epochs. 


The  Third  Epoch 


For  the  third  epoch  we  used  plates  exposed  by  Richard  Smart,  again  on  the  26-inch 
refractor  of  the  Leander  McCormick  Observatory.  For  this  epoch  we  recorded  each 
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Figure  10:  The  Leandcr  McCormick  Refractor 
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chosen  region  only  an  average  of  3 times  on  a single  plate  as  opposed  to  4.  However 
we  covered  97  regions,  allowing  approximately  the  same  coverage  but  with  greater 
overlap.  Figure  11  illustrates  the  orientation  of  the  regions  chosen. 

Both  the  second  and  third  epoch  were  observed  using  a large,  2.5  magnitude 
objective  grating.  This  meant  that  most  bright  ( i.e.,  brighter  than  10th  magnitude) 
stars  formed  up  to  5 measurable  stellar  images.  Later  on  we  will  discuss  how  these 
images  were  used  to  control  the  possible  magnitude  effects.  Table  4 summarizes  the 
region  data  listed  in  the  same  format  as  the  second  epoch  data. 
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Table  4;  Third  Epoch  McCormick  Observatory  Plates  Data 


MC# 

Ooh 

o 

o 

Epoch 

EST 

HA 

s 

T 

Temperature 

141225 

5.28 

-1.3 

1991.9 

12:45 

0:25W 

2 

2 

56 

141227 

5.26 

-1.2 

1991.9 

13:14 

0:10W 

3 

2 

57 

141229 

5.38 

-1.1 

1991.9 

14:04 

0:25E 

2 

2 

57 

141231 

5.65 

-1.7 

1991.9 

14:36 

0:50E 

2 

2 

57 

141247 

5.25 

-1.5 

1991.9 

12:31 

0:50E 

3 

5 

55 

141249 

5.38 

-1.5 

1991.9 

13:00 

0:30E 

3 

4 

55 

141251 

5.46 

-1.5 

1991.9 

13:38 

0 : 05W 

3 

4 

55 

141253 

5.52 

-1.5 

1991.9 

14:05 

0 : 30W 

3 

3 

55 

141263 

5.35 

-0.2 

1991.9 

12:34 

0 : 40W 

3 

5 

33 

141265 

5.41 

-0.2 

1991.9 

13:03 

0.-20W 

3 

5 

33 

141267 

5.48 

-0.2 

1991.9 

13:34 

0:08E 

4 

5 

30 

141269 

5.55 

-0.2 

1991.9 

13:56 

0:36E 

4 

5 

31 

141279 

5.28 

-0.6 

1991.9 

12:51 

0:20E 

4 

5 

33 

141281 

5.35 

-0.6 

1991.9 

13:18 

0:04W 

4 

5 

33 

141283 

5.41 

-0.6 

1991.9 

13:42 

0:23W 

4 

5 

33 

141285 

5.48 

-0.6 

1991.9 

14:06 

0 : 4 5W 

4 

5 

32 

141341 

5.36 

0.0 

1991.9 

11:59 

0:06E 

3 

5 

42 

141342 

5.42 

0.0 

1991.9 

12:17 

0 : low 

3 

5 

42 

141343 

5.47 

0.0 

1991.9 

12:32 

0:21W 

3 

5 

42 

141344 

5.53 

0.0 

1991.9 

13:00 

0:45W 

2 

5 

41 

141345 

5.25 

-1.8 

1991.9 

11:10 

0:44E 

2 

4 

37 

141346 

5.31 

-1.8 

1991.9 

11:37 

0:20E 

2 

4 

37 

141347 

5.36 

-1.8 

1991.9 

11:57 

0:03E 

2 

4 

37 

141348 

5.42 

-1.7 

1991.9 

12:13 

0 : low 

2 

4 

36 

141349 

5.48 

-1.8 

1991.9 

12:36 

0:28W 

2 

4 

36 

141350 

5.54 

-1.8 

1991.9 

12:46 

0 : 36W 

2 

4 

36 

141351 

5.60 

-1.7 

1991.9 

13:04 

0:50W 

2 

4 

36 

141352 

5.72 

-1.8 

1991.9 

13:20 

0:58W 

2 

4 

36 

141353 

5.34 

-2.2 

1991.9 

11:25 

0:30E 

4 

4 

32 

141354 

5.39 

-2.1 

1991.9 

11:42 

0:15E 

4 

4 

32 

141355 

5.45 

-2.1 

1991.9 

11:59 

0:03E 

4 

4 

32 

141356 

5.51 

-2.1 

1991.9 

12:24 

0 : 16W 

4 

4 

32 

141357 

5.57 

-2.1 

1991.9 

12:47 

0 : 34W 

4 

4 

32 

141358 

5.63 

-2.1 

1991.9 

12:57 

0 : 46W 

3 

2 

32 

141368 

5.26 

-2.4 

1992.0 

11:02 

0 : 40E 

2 

4 

30 

141369 

5.31 

-2.4 

1992.0 

11:29 

0:16E 

2 

4 

30 

141370 

5.36 

-2.4 

1992.0 

11:55 

0 : 06W 

2 

4 

30 

141371 

5.42 

-2.4 

1992.0 

12:19 

0:28W 

2 

3 

29 

141372 

5.48 

-2.4 

1992.0 

12:36 

0 : 40W 

2 

3 

29 

141373 

5.54 

-2.4 

1992.0 

12:57 

0 : 58W 

2 

3 

29 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 

61 

62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 
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MC# 

Ooh 

Epoch 

EST 

HA 

s 

T 

Temperature 

141414 

5.28 

-2.7 

1992.0 

09:17 

0:37E 

3 

3 

32 

141415 

5.35 

-2.7 

1992.0 

09:34 

0 : 23E 

3 

3 

32 

141416 

5.40 

-2.7 

1992.0 

10:20 

0 :20W 

2 

3 

31 

141417 

5.46 

-2.7 

1992.0 

10:40 

0:36W 

2 

3 

31 

141418 

5.69 

-2.7 

1992.0 

11:02 

0 : 44w 

2 

3 

31 

141419 

5.30 

-3.0 

1992.0 

08:55 

0:56E 

4 

3 

25 

141420 

5.37 

-3.0 

1992.0 

09:20 

0:  34E 

4 

3 

25 

141421 

5.42 

-3.0 

1992.0 

09:39 

0:20E 

4 

3 

25 

141422 

5.49 

-3.0 

1992.0 

09:59 

0:03E 

4 

3 

25 

141423 

5.54 

-3.0 

1992.0 

10:17 

0 : IIW 

3 

3 

22 

141424 

5.61 

-3.0 

1992.0 

10:34 

0:25W 

3 

3 

22 

141425 

5.66 

-3.0 

1992.0 

10:49 

0 : 37W 

3 

3 

22 

141426 

5.72 

-3.0 

1992.0 

11:10 

0 : 55W 

3 

3 

22 

141429 

5.35 

-1.0 

1992.0 

08:54 

0:47E 

3 

3 

22 

141430 

5.40 

-1.0 

1992.0 

09:09 

0 : 35E 

3 

3 

22 

141431 

5.28 

-1.0 

1992.0 

09:30 

0:09E 

4 

3 

28 

141432 

5.45 

-0.9 

1992.0 

09:48 

0:00E 

4 

3 

28 

141433 

5.52 

-0.9 

1992.0 

10:09 

0 : 17W 

4 

3 

28 

141434 

5.57 

-0.9 

1992.0 

10:26 

0 : 31W 

4 

3 

28 

141435 

5.62 

-0.9 

1992.0 

10:35 

0-.36W 

4 

3 

28 

141444 

5.34 

-1.5 

1992.1 

08:21 

0:53E 

3 

4 

34 

141445 

5.36 

0.6 

1992.1 

08:38 

0:37E 

3 

4 

34 

141446 

5.42 

0.6 

1992.1 

08:57 

0:22E 

3 

4 

34 

141447 

5.48 

0.6 

1992.1 

09:21 

0 : 02W 

3 

4 

32 

141448 

5.45 

-0.3 

1992.1 

09:36 

0 : 16W 

3 

4 

32 

141449 

5.51 

-2.7 

1992.1 

10:00 

0:37W 

3 

4 

32 

141450 

5.57 

-2.7 

1992.1 

10:15 

0:47W 

3 

4 

32 

141455 

5.34 

-3.3 

1992.1 

08:30 

0:36E 

3 

4 

42 

141456 

5.40 

-3.3 

1992.1 

08:46 

0:23E 

3 

4 

42 

141457 

5.45 

-3.3 

1992.1 

09:03 

0:10E 

3 

4 

41 

141458 

5.51 

-3.3 

1992.1 

09:24 

0 : 09w 

3 

4 

41 

141459 

5.57 

-3.3 

1992.1 

09:38 

0 : 18W 

3 

4 

41 

141460 

5.63 

-3.3 

1992.1 

09:55 

0:32W 

3 

4 

41 

141461 

5.60 

-2.4 

1992.1 

10:15 

0:54W 

3 

4 

41 

141471 

5.48 

-1.2 

1992.1 

09:20 

0 : 09W 

3 

3 

45 

141472 

5.54 

-1.2 

1992.1 

09:34 

0:20W 

3 

3 

45 

141473 

5.60 

-1.2 

1992.1 

09:45 

0:27W 

3 

3 

45 

141474 

5.57 

-1.5 

1992.1 

10:00 

0:45W 

3 

3 

45 

141475 

5.63 

-1.5 

1992.1 

10:14 

0:55W 

3 

3 

45 

141490 

5.39 

0.9 

1992.1 

07:33 

0:06E 

4 

4 

44 
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# 

MC# 

Ooh 

o 

o 

Epoch 

EST 

HA 

s 

T 

Temperature 

81 

141491 

5.45 

0.9 

1992.1 

07:48 

0 : 06W 

4 

4 

44 

82 

141492 

5.54 

-0.6 

1992.1 

08:07 

0:20W 

4 

4 

44 

83 

141493 

5.60 

-0.6 

1992.1 

08:26 

0:35W 

4 

4 

44 

84 

141494 

5.63 

-2.7 

1992.1 

08:38 

0:46W 

4 

4 

44 

85 

141511 

5.34 

-3.9 

1992.2 

07:12 

0 : 44w 

4 

4 

62 

86 

141512 

5.40 

-3.9 

1992.2 

07:26 

0 : 55w 

4 

4 

62 

87 

141513 

5.45 

-3.9 

1992.2 

07:35 

0 : 59w 

4 

4 

62 

88 

141514 

5.63 

-3.9 

1992.2 

07:48 

1:02W 

4 

4 

62 

89 

141515 

5.72 

-2.4 

1992.2 

08:00 

1:07E 

4 

4 

62 

90 

141522 

5.37 

-3.6 

1992.2 

07:17 

0:50W 

2 

4 

34 

91 

141523 

5.44 

-3.6 

1992.2 

07:25 

1 : OOW 

2 

4 

34 

92 

141524 

5.64 

-3.6 

1992.2 

07:37 

0:55W 

2 

4 

34 

93 

141525 

5.68 

-2.1 

1992.2 

07:42 

1:05W 

2 

4 

34 

Measuring  the  Plates 

We  measured  the  plates  on  the  University  of  Minnesota  Automated  Plate  Scanner 
(APS)  at  Minneapolis,  which  was  originally  designed  as  a proper  motion  measuring 
engine  for  Luyten’s  automated  proper  motion  surveys.  It  has  been  refitted  with  modem 
electronics  and  is  currently  being  used  for  digitizing  the  Palomar  Sky  Surveys. 

The  machine  is  shown  in  two  pictures  in  figures  12  and  13  (figure  12  shows  the 
overall  machine.)  It  can  measure  two  plates  at  any  one  time  on  the  two  measuring 
levels.  For  measuring  the  second  and  third  epoch  plates  we  used  just  the  upper  level. 
Figure  13  shows  the  laser  beam  that  measures  the  plate  density  as  the  whole  housing 
scans  across  the  plate.  Inside  this  housing  is  an  eight  sided  prism  that  rotates,  with 
the  speed  of  the  rotation  setting  the  scanning  speed.  It  scans  out  12mm  strips  with  12 
micron  diameter  spots  in  y and  simultaneously  the  strip  is  moved  in  x.  Thus  knowing 
the  speed  of  rotation,  the  motion  in  x,  and  the  density  at  any  point  a digitized  map  of 
the  plate  can  be  made. 
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Figure  12:  The  Automated  Plate  Scanner(Pennington  et  al,  1993) 
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Figure  13:  The  APS  Scanning  Laser  {Pennington  et  ai,  1993) 
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Siam 

Figure  14:  APS  Isodensi  tome  trie  scan  (Pennington  et  ai,  1993) 

In  addition  to  differing  speeds  (and  as  such  differing  resolutions)  there  are  three 
different  kinds  of  scanning  modes  for  the  APS:  isodensitometric,  densitometric  and  a 
hybrid  of  these  two  modes.  For  measuring  these  plates  we  used  the  isodensitometric 
mode.  In  this  mode  the  diffraction  pattern  of  the  laser  beam  as  it  scans  is  focussed 
on  a photodetector.  Every  time  the  image  attains  either  a preset  or  a dynamically  set 
threshold  density,  the  position  is  read  to  a precision  of  0.366  microns.  Therefore  each 
stellar  image  will  have  two  points  for  each  scan  across  it,  the  ingress  and  egress  of 
the  image.  Figure  14  shows  a typical  stellar  image  scanned  in  isodensitometric  mode 
with  four  threshold  levels. 

An  ellipse  is  fitted  to  this  isodensitometric  outline.  This  best  fit  ellipse  then  defines 
the  image  and  its  center  defines  the  position  of  the  image.  The  final  output  consists 
of  six  parameters  for  each  image:  x and  y coordinates  of  the  ellipse  center,  diameter, 
ellipticity,  theta  angle  (angle  of  the  best  fit  ellipse’s  semimajor  axis  to  the  x axis),  and 
a value  that  measured  the  goodness  of  fit  of  the  ellipse.  For  a typical  plate  there  may 
be  5000  images  of  which  only  1000  were  real  stellar  images  or  grating  images. 

This  machine  is  extremely  fast  and  allowed  us  to  measure  all  the  plates  twice  in 
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the  span  of  two  weeks.  For  the  first  scan  we  only  used  one  level  threshold;  for  the 
second  scan  we  used  two  thresholds.  So  for  each  plate  we  have  three  estimates  of  the 
images  on  that  plate. 

Each  adjacent  stripe  has  a 1mm  overlap  region.  Images  in  this  overlap  region 
that  have  matching  parameters  are  assumed  to  be  the  same  image  and  then  used  to 
tie  stripes  together.  This  is  because  the  machine  would  tend  to  drift  by  perhaps  2-3 
microns  between  stripes  and  the  use  of  an  overlap  allowed  a continual  mapping  of  the 
measuring  frame  to  a fixed  frame. 

Figure  15  shows  the  predicted  stellar  positions  for  the  regions  141368,  141369  and 
141370.  Each  region  was  exposed  three  times  on  the  same  emulsion.  One  emulsion 
was  used  due  to  their  high  cost,  it  also  made  the  scanning  easier  as  we  scanned  the 
whole  plate  in  one  ‘sweep’.  Each  plate  is  filed  under  the  McCormick  number  of  the 
first  exposure. 

Each  star  with  a cross  is  a reference  star;  they  are  brighter  on  average  than  the 
field  stars.  The  grating  images  of  the  brighter  stars  (m<10)  are  used  to  ensure  that 
any  magnitude  effects  are  equally  modelled  for  the  whole  of  the  stellar  magnitude 
range.  With  ideal  magnitude-independent  imaging,  the  primary  stellar  image  will  be 
at  the  average  of  the  first  and  second  order  images.  Coma  and  magnitude  effects  will 
cause  deviations  from  an  exact  agreement  of  the  average  of  the  grating  images  and 
the  primary  images.  By  examining  this  deviation  we  can  very  accurately  model  the 
magnitude  dependent  terms.  (There  is  no  better  way  to  do  it.) 

There  are  different  ways  of  examining  this  deviation.  In  this  analysis  for  every 
measured  objective  grating  pair  we  assigned  a fictitious  star  to  their  average.  Figure  16 


141368x70 
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Figure  16:  Predicted  Positions  for  Plate  141368x70 
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shows  an  example  of  a difference  that  might  be  expected.  The  X marks  the  average 
of  the  two  first  order  grating  images.  At  this  point  we  ‘create’  a fictitious  star  that 
has  all  the  same  parameters  as  the  primary  star  except  its  magnitude  is  reduced  by  the 
grating  constant.  This  is  then  entered  into  the  overlap  as  an  extra  observation  and  the 
information  contained  in  its  position  and  magnitude  will  be  enforced  by  the  overlap 
least-squares  reduction. 


Grating  images 


image. 


Figure  16:  The  difference  between  a Grating  Image  Average  and  the  Primary  Image  Center 


This  method  may  increase  the  variance  of  the  solution,  but  we  feel  the  possible 
increase  in  variance  due  to  extra  observations  will  be  outweighed  by  the  potentially 
stronger  overlap  and  the  ability  to  solve  for  all  the  parameters  simultaneously. 

Figure  17  shows  the  scan  of  plate  141368.  This  shows  one  drawback  of  using  a 
diffraction  grating.  The  large  ellipses  are  blended  diffraction  images  that  the  software 
has  been  unable  to  differentiate.  These  images  were  not  used  in  the  stellar  positions 
and  this  results  in  a minor  loss  of  data.  As  we  have  3-4  exposures  of  each  region,  and 
because  a particular  star  may  appear  in  perhaps  4 regions  then  the  loss  of  1 or  2 of  its 
measurements  will  not  greatly  affect  the  final  positions. 

To  find  the  location  of  the  bright  stars  that  are  differentially  subject  to  image 
blending  in  the  second  set  of  scans  (the  reverse  scans),  we  used  two  thresholds.  Figure 
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18  shows  the  same  plate  scanned  again  in  the  upper  threshold.  As  the  figure  shows 
we  have  lost  the  dimmer  stars  and  grating  images,  but  the  brighter  stars  have  been 
resolved.  This  was  treated  as  a separate  exposure  of  the  same  region,  again  allowing 
a strong  overlap  and  actually  allowing  an  independent  check  of  the  plate  parameters; 
they  should  be  the  same  for  the  second  scan. 

Summary  of  Observational  Material  and  Measured  Images 

To  summarize,  the  observational  material  basically  consists  of  photographic  plates 
exposed  at  three  epochs.  Below  are  the  observation  statistics  for  all  three  epochs: 


First  Epoch  Plate  Information 

Observatory  and  epoch 
Number  of  Regions 
Number  of  Exposures 
Plate  scale  of  telescope 
Average  Positional  precision 

San  Fernando  and  Algiers,  1900 

31 

31 

607mm 

4/fm 

Second  Epoch  Plate  Information 

Observatory  and  epoch 
Number  of  Regions 
Number  of  Exposures 
Plate  scale  of  telescope 
Average  Positional  precision 

McCormick  Observatory,  1956 
67 
261 

227mm 

2;/m 

Third  Epoch  Plate  Information 

Observatory  and  epoch 
Number  of  Regions 
Number  of  Exposures 
Plate  scale  of  telescope 
Average  Positional  precision 

McCormick  Observatory,  1992 

93 

293 

227mm 

2//m 

R141368, 
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Figure  17:  Scan  of  Plate  141368x70  Low  Threshold 


R141368 
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Figure  18:  Scan  of  Plate  141368x70  High  Threshold 


CHAPTER  5 
FINDING  A MODEL 


In  the  discussion  on  Finding  Stellar  Coordinates  we  used  as  an  example  the  six- 
constant  model.  However,  as  we  also  discussed,  to  account  accurately  for  the  complex 
relationship  between  the  observations  (x,y)  and  the  standard  coordinates  (^,7/)  a greater 
number  of  terms  is  needed.  The  model  will  invariably  not  account  for  all  contributing 
effects.  There  comes  a point  when  model  improvement  will  require  more  effort  than 
the  potential  improvement  in  accuracy  warrants.  Also,  applying  a complicated  model 
known  to  be  correct  but  which  requires  the  estimation  of  a relatively  large  number 
of  parameters  may  depress  the  ultimate  accuracy  due  to  generating  a large  parameter 
variance,  (cf.  Eichhorn  and  Williams,  1963). 

For  the  choice  of  a model  there  is  no  unambiguous  technique  that  will  give  black 
and  white  indications  of  which  parameters  to  include.  The  procedure  is  more  of  an 
art  than  a mathematical  reduction.  In  addition  to  the  numerical  tests  of  significance  an 
understanding  of  the  physical  situation  aids  the  analysis  significantly. 

The  discussion  will  be  split  into  sections:  previous  studies,  position  residuals 
analysis,  parameter  variation  analysis,  and  solution  variances.  There  will  also  be  a 
separate  discussion  of  single  and  overlapping  plate  solutions,  this  is  important  because 
the  noise  in  the  single  plate  solution  may  hide  some  effects  and  the  stronger  tie  in  the 
overlap  may  also  hide  some  different  effects.  These  sections  are  chosen  to  simplify 
the  discussion.  The  final  choice  was  based  on  an  examination  all  of  these  indicators 
simultaneously,  while  trying  out  different  models. 
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Previous  Work 


The  original  publications  of  the  Astrographic  Catalogue  included  the  basic  six- 
constant  models.  This  was  because  the  accuracy  required  only  warranted  a six-constant 
model  and  the  number  of  reference  stars  per-square  degree  never  really  allowed  a larger 
model.  Using  the  ACRS  there  is  an  average  of  32  reference  stars  per  plate.  This 
will  give  us  only  a parameter- to-observation  ratio  of  two  in  a 16-constant  single  plate 
solution.  This  ratio  is  much  higher  in  the  overlap  where  there  are  on  average  200 
images  per  plate,  many  of  which  will  act  as  observations.  Therefore  for  reasonable  size 
models  by  using  an  overlap  solution  we  are  no  longer  limited  to  a small  model  and  we 
can  examine  what  other  possible  terms  may  be  included. 

Previous  work  on  zones  of  the  AC  have  used  simple  six-constant  models  to  ten- 
constant  models  (Eichhom  and  Gatewood,  1967)  all  the  way  up  to  a combination  of 
the  six-constant  model  with  twenty  constrained  constants  (Gunther  and  Kox,  1971). 
Various  investigations  have  found  that  AC  material  contains  evidence  for  magnitude 
terms,  color  terms,  coma  terms  and,  radial  terms. 

In  this  analysis  we  deal  with  the  Algiers  and  the  San  Fernando  zone.  During 
the  measuring  of  the  Algiers  zone  no  magnitude  term  was  found  (Eichhom  and  Clary, 
1974).  In  the  San  Fernando  zone  the  measurements  showed  no  signs  of  radial  distortion, 
aberration  or  refraction.  (It  is  worth  noting  however  that  high  discrepancies,  up  to 
25//,  between  rectangular  coordinates  of  the  same  image  were  tolerated,  and  thus  these 
effects  could  easily  be  lost  in  the  noise.) 

The  astrometric  properties  of  the  McCormick  Refractor  were  examined  by  Rus- 
sell (1976).  She  examined  data  from  15  large  astrometric  telescopes  for  evidence  of 
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model  terms  in  magnitude,  coma,  magnitude  squared,  color,  color  magnification  (color- 
coordinate  product),  tilt,  refraction  (position  cross  product),  distortion  and,  magnitude- 
distortion  product.  The  McCormick  Telescope  showed  significant  terms  in  the  magni- 
tude, coma,  tilt,  refraction,  and  the  magnitude-distortion  term.  When  we  finalize  our 
model  choice  we  will  refer  to  the  previous  work  to  see  if  they  agree  with  previous 
investigations. 

Positional  Residuals 

In  this  section  we  will  analyze  the  output  from  six-constant  solutions.  The  objective 
is  to  use  a very  basic  model  and  see  what  extra  terms  the  position  residuals  and  errors 
are  correlated  with.  We  analyze  21  extra  terms  by  graphing  the  residuals  of  our  solution 
vs  their  value.  The  21  extra  terms  analyzed  were 

{m-  {m)),(m-  {m)f,{m-{m))x,{7n-  (m})y, 

(m  - {m))x~,  (m  - (m))y-,  (m  - (m)fx,  (m  - {m)fy, 

(m-  (?n))(x-  + y‘),(?n-  (m))x(x^  + y^),  (5.1) 

(m-  (m))?/(a;2  + 2/'-),(,7i-  {in))(xy),xy,x^,y-, 
x^y,  xy-,  (x^  -f  7J-)  ,x{x-  + y^),y{x'^  + y-),  (x^  -|-  y^)^ . 

In  the  analysis  of  each  term  we  examined  a host  of  various  graphs.  It  would  be 
counter  productive  to  reproduce  all  of  the  graphs,  examples  are  reproduced  for  each 
technique. 
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Single  Plate  Residuals 

Consider  the  results  from  a single  plate  solution.  For  the  reference  stars  the  residuals 
can  be  found  and  are  of  the  form 


where  the  x,  y,  ( and  rj  are  the  inputted  values  and  the  parameter  matrix  is  the  result 
of  the  reduction.  These  residuals  can  then  be  examined  to  see  which  other  terms  in 
the  stellar  parameters  may  be  significant.  In  figures  19  and  20  we  graph  the  residuals 
(Ax  and  Ay)  of  a six-constant  single  plate  solutions  vs  the  different  extra  terms  listed 
above.  The  residuals  are  for  27  plates  of  the  1900  epoch  material;  all  plates  with  over 
25  reference  stars  were  included. 

The  residual  on  the  y axis  in  each  graph  is  limited  to  three  standard  deviations,  the 
dotted  line  represents  two  standard  deviations.  There  is  a large  amount  of  scatter  in 
the  points.  This  is  expected  as  ideally  they  would  be  randomly  distributed  reflecting 
the  problems  in  the  reference  data  and  measuring  process.  There  may  also  be  some 
misidentifications  in  the  sample  as  well  which  will  add  to  the  noise. 

A 3rd  order  polynomial  was  fitted  to  the  points,  for  the  purpose  of  examining 
trends  in  the  curve.  Graphs  that  were  similar,  both  in  scale  and  trend,  were  noted.  For 
example,  if  we  found  a trend  in  a term  that  was  similar  in  both  the  x and  y residual  then 
we  can  assume  the  parameter  for  this  term,  should  we  include  it,  would  be  the  same. 

From  figure  19  we  can  make  the  following  conclusions:  the  residuals  are  correlated 
with  both  the  m (magnitude)  and  m^  quantities,  the  correlation  is  stronger  and  consistent 
in  the  Ay,  and  not  as  strong  or  consistent  in  the  Ax.  There  is  no  strong  evidence  for 


(5.2) 
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Figure  19:  Six-Constant  Single  Plate  Residuals 
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Figure  20:  Six-Constant  Single  Plate  Residuals 
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a correlation  in  the  Ax  vs  mx,  my  or  mx^.  In  the  Ay  vs  mx  and  my  there  is  evidence 
for  a correlation  and  only  slight  evidence  in  the  Ay  vs  mx^. 

From  figure  20  we  can  see  that  the  Ax  is  not  consistently  correlated  in  any  of  the 
terms,  though  it  may  be  in  the  mx(x^+y^)  term.  The  Ay  shows  a strong  correlation  in 
my^,  xm^,  and  ym^  and  a weaker  correlation  in  the  mx(x^+y^)  and  my(x^+y^)  terms. 

These  correlations  imply  that  the  plates  suffered  from  a guiding  error  and  some 
coma,  mainly  in  the  y coordinate.  Based  on  these  graphs,  for  the  1900  model,  we 
should  include  terms  in  the  x coordinate  that  fit  m,  m^  and  prehaps  mx,  and  in  the  y 
coordinate  terms  that  fit  m,  m^,  mx,  my,  xm^,  and  ym^.  In  the  actual  development  we 
would  not  do  this  straight  away.  First  we  would  examine  the  corresponding  six-constant 
model  overlap  residuals  (as  in  the  next  section),  and  then  if  there  were  still  the  same 
correlations  we  would  fit  a magnitude  term  and  one  or  two  of  the  other  terms.  Then 
we  reexamine  the  remaining  terms  and  add  terms  as  needed. 

Overlap  Residuals 

In  the  single  plate  method  because  each  solution  was  only  dependent  on  20-30 
reference  stars,  the  solutions  could  be  loose  and  therefore  some  of  the  trends  could 
have  been  hidden  in  the  noise.  In  the  overlapping  plate  method  the  solution  is  much 
stronger  because  of  the  restriction  that  a star  have  the  same  position  from  plate  to  plate. 
Therefore,  even  though  we  are  only  using  a six-plate  model  some  effects,  such  as  radial 
ones,  will  be  accounted  for  by  the  overlap. 

For  example,  the  images  of  a given  star  at  the  centers  of  one  plate  must  have  the 
same  position  as  an  image  of  the  same  star  that  may  be  in  the  corner  on  another  plate 
at  any  give  epoch.  Any  radial  effects  will  be  canceled  out  when  we  find  the  mean 
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weighted  position.  When  we  find  the  residual  of  the  position  this  effect  is  slightly 
nullified  because  we  are  returning  again  to  single  images  rather  than  mean  weighted 
average  star  positions. 

Now  we  examine  the  residuals  from  a six-constant  overlap  as  we  did  in  the  single 
plate  example  to  see  if  there  are  any  improvements  or  any  other  problems  highlighted. 
Figures  21  and  22  show  examples  of  residual  graphs  for  a six-constant  overlap  of  the 
1900  material. 

In  figure  21  we  again  see  the  same  positive  correlation  in  the  m and  m^  terms, 
though  it  is  damped  slightly.  There  does  however  seem  to  be  no  correlation  the  other 
terms  except  for  the  my  term  which  shows  a strong  parabolic  shape;  this  implies  a 
very  real  trend. 

In  figure  22  we  see  some  trends  that  we  had  thought  were  not  consistent,  like  that 
in  the  Ax  vs  mx(x^-i-y^)  and  my(x^-f-y^)  terms  and  we  also  see  a new  similar  feature  in 
the  xm^  term.  If  we  decide  to  fit  the  mx(x^-t-y^)  and  my(x^-i-y^)  for  either  Ax  or  Ay  we 
can  keep  the  same  parameters  for  them  as  the  correlation  is  identical.  The  only  other 
correlation  is  in  the  Ay  vs  ym^  term  which  repeats  the  correlation  seen  in  the  single 
plate  analysis  and  should  be  modelled. 

Summary 

From  this  analysis  we  can  highlight  those  terms  we  feel  important  for  further 
analysis.  After  the  more  sophisticated  model  is  run,  another  examination  of  the  residuals 
may  highlight  more  correlations  and  these  will  have  to  be  included.  An  analysis  of  the 
parameter  variation  in  the  next  section  shows  how  we  can  highlight  parameters  that 
may  be  correlated  but  are  not  significant  to  the  reduction  and  therefore  can  be  removed. 
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Figure  21:  Six-Constant  Overlap  Residuals 
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Figure  22:  Six-Constant  Overlap  Residuals 
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From  our  residual  analysis  we  found  that  all  three  epochs  showed  significant  first 
order  magnitude  effects  and  so  the  minimum  model  would  be  the  eight  constant  model: 

_ /'^TjlOOOm  — (m)  0 \ 

\000^7?1  0 m—{m))’ 

Now  for  each  of  the  epochs  we  found  different  effects  and  they  are  discussed  separately. 

Epoch  1900:  An  analysis  of  the  residuals  showed  the  terms  were  correlated  mainly 
with  the  Ay  residual.  This  Ay  varied  so  strongly  in  one  term  the  Ax  residual  is  not 
modelled  but  the  Ay  is.  A very  strong  correlation  was  also  found  in  the  first  and  second 
order  magnitude  terms  that  could  possibly  be  modeled  by  only  two  parameters.  There 
is  a consistent  and  similar  variation  in  the  mx(x^+y^)  and  my(x^+y^)  terms  indicating 
they  should  only  be  modeled  in  y and  that  it  may  be  the  m(x^+y^)  term  that  needs 
modelling.  We  found  a similar  but  opposite  (e.g.  acting  on  the  opposite  residual)  in  the 
xy  and  mxy  terms,  signifying  these  can  be  modeled  by  only  two  parameters,  similarly 
in  the  xy^  and  x(x^+y^)  terms.  Finally  correlations  were  different  in  the  Ay  and  Ax 
vs.  (x^+y^),  hence  two  parameters  were  needed  for  this  one  radial  term. 

After  residual  analysis  we  decided  to  try  the  eight  constant  model  plus  the  following 
eight  terms: 


a 

b' 

c' 

e 


(5.3) 


(5.4) 
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Epoch  1956:  For  the  1956  material  there  was  very  little  evidence  for  magnitude 
and  tilt  terms.  The  coma  and  magnitude  distortion  terms  seem  to  dominate,  but  even 
these  cause  very  little  variation.  After  the  residual  analysis  we  decided  to  examine  the 
following  terms: 

(m  - 0 (m  - {m))7(($^ + rj^)  0 

0 (m  - (m))^ 0 (m  - 


(5.5) 


{m-(m)f^  0 {m  - {m)f  + T)^) 

0 {m-{m))^r}  (m  - (m))^ 

Epoch  1992:  Again  dominated  by  coma  and  magnitude  terms;  in  this  case  however 
the  graphs  are  similar  and  we  have  modelled  these  terms  with  single  parameters.  There 
was  some  slight  slope  in  the  residual  plots  for  the  tilt  terms  and  origin-distance  squared 
terms  but  not  very  pronounced.  After  residual  analysis  we  decided  to  examine  the 
following  terms: 

(m  - (m))^(^2  ^ ^ ^2^ 

{m  - + r]'^)  (m  - (m))^ ^ 


(5.6) 


0 (m-(?n))^/  (m  - (m))^(^2  ^ j^2) 

Note  that  even  though  the  residuals  may  not  imply  there  is  a magnitude  term  we 
have  included  one.  This  is  because  the  magnitude  effect  has  a strong  tie  to  the  particular 
guiding  in  an  exposure;  all  the  plates  were  hand  guided,  and  therefore  the  term  is  likely 
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to  be  randomized  over  the  entire  plate  system.  As  such  it  may  not  show  up  in  this 
kind  of  analysis;  on  the  other  hand  if  it  is  negligible  then  an  analysis  of  the  parameter 
variation  will  highlight  this. 

Two  things  to  note  about  the  models:  the  astrographic  telescopes  do  not  show 
particular  coma  or  magnitude  distortion  terms,  but  do  show  some  radial  distortion  and 
tilt  terms.  The  1956  and  1992  models  agree  quite  well,  as  one  would  expect  when  using 
the  same  telescope.  They  also  contains  terms  in  magnitude  distortion  and  coma  found 
by  Russell  (1976),  though  there  was  no  strong  evidence  for  tilt  terms  as  she  found. 

Parameter  Variation 

Now  we  have  picked  the  terms  we  feel  are  correlated  with  the  residuals  and  should 
therefore  be  modeled.  We  should  now  check  to  make  sure  all  the  terms  are  needed.  It 
would  be  quite  possible  for  a term  to  be  highly  correlated  but  still  have  a negligible 
effect.  Removing  extra  terms  is  important  for  two  considerations: 

1.  As  one  increases  the  number  of  parameters  their  variance  will  also  increase.  This 
will  propagate  an  error  directly  into  the  calculated  positions.  Eichhorn  and  Williams 
showed  that  parameter  variance  can  actually  make  a correct  model  give  less  accuracy 
than  a smaller  approximate  model. 

2.  In  the  physical  process  of  computation  every  calculation  results  in  a small  round- 
off error.  If  we  consider  a overlap  of  30  plates  with  12  parameters,  then  the 
matrix  of  nonnal  equations  has  dimension  (360  x 360).  If  we  now  include  2 more 
parameters  the  dimension  increases  to  (420  x 420).  A good  rule  of  thumb  is  that 
the  inversion  of  a matrix  will  require  multiplications  (where  N is  the  order  of  the 
matrix).  Therefore,  continuing  with  our  example  the  inclusion  of  2 extra  parameters 
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will  require  (420^  - 360^)  = 27,432,000  extra  calculations,  nearly  half  as  many 
multiplications  as  those  needed  for  the  12  constant  model.  The  loss  of  accuracy 
incurred  by  round  off  error  in  the  increased  computations  may  well  outweigh  the 
potential  accuracy  increase  due  to  using  a better  model. 

The  first  examination  must  be  made  with  the  single  plate  analysis,  because  including 
too  many  terms  in  the  overlap  will  cause  the  inversion  process  to  ‘blow  up’  and 
the  matrix  will  not  be  invertible.  This  examination  quantifies  the  graphically  based 
subjective  choices  made  above.  It  will  enable  a simplifying  of  the  potential  models  by 
highlighting  parameters  that  have  negligible  effects  (even  highly  correlated  ones)  and 
parameters  that  do  not  change  from  plate  to  plate. 

The  parameter’s  mean  and  standard  deviation  will  indicate  how  well  determined 
and  significant  the  parameter  is.  The  mean  parameter  and  the  largest  ^ and  t]  values 
will  provide  an  upper  bound  on  the  term’s  contribution  . For  the  AC  the  largest 
standard  coordinate  value  is  if,r/=0.038  radians  and  for  the  McCormick  data  ^=0.018 
and  r/=0.013radians. 

The  errors  of  the  individual  parameters  will  also  be  important  to  examine.  These 
will  be  given  by  the  square  root  of  the  diagonal  term  of  the  covariance  matrix.  In  a 
system  of  equations  represented  by  Aa  = x„  the  covariance  matrix,  C,  is  proportional 
to  the  inverse  of  the  matrix  A.  Referring  back  to  the  single  plate  equation  (3.23  ) 


here  the  matrix 


is  proportional  to  the  covariance  matrix. 
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Finally,  terms  should  be  examined  to  see  if  they  are  correlated.  To  do  this  we 
calculate  the  correlation  parameter:  where  Cv^^  represents  the  term  of  the 

covariance  matrix. 

Table  5 summarizes  all  these  quantities  for  the  single  plate  solutions  using  the 
models  outlined  in  5.4,  5.5,  and  5.6.  As  in  the  single  plate  solution  we  only  use  the 
reference  stars,  these  parameters  will  not  be  very  well  determined,  but  their  general 
variation  will  provide  information  for  the  removal  or  constraining  of  parameters  as 
required.  Listed  are 


#: 

Parameter  number,  7 and  8 represent  first  order  magnitude  terms,  the 
other  parameters  are  as  given  in  equations  5.4,  5.5,  and  5.6. 

Mean: 

Mean  parameter  value,  if  this  is  very  large  it  means  the  parameter 
will  be  dominating  the  calculation  of  the  whole  parameter  vector. 

a\ 

Standard  deviation  of  the  parameter  value,  if  this  is  greater  than  the 
mean  value  then  one  must  consider  if  the  parameter  can  be  ignored. 

Range: 

Range  of  the  parameter,  indicates  if  there  are  some  aberrant  values. 

Mean  e: 

Mean  parameter  error,  this  is  the  square  root  of  the  covariance  matrix 
diagonal  term  averaged  for  all  the  plates. 

ot- 

Standard  deviation  of  the  error. 

Mean  t ! a\ 

Ratio  of  the  mean  error  to  the  parameter  standard  deviation.  If  this 
is  equal  to  one  then  the  parameter  is  varying  within  its  error. 

Mean  / Mean  e: 

Ratio  of  the  mean  parameter  value  to  the  mean  error.  If  this  is  equal 

to  one  then  the  parameter  is  within  its  error,  zero. 
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a tier.  Ratio  of  error  to  parameter  standard  deviation.  If  this  is  one,  65%  of 

the  values  are  within  their  error  the  mean  value. 

Term  Mean:  The  highest  value  of  the  model  term  if  you  use  the  mean  value  of 

the  parameter. 

Term  Range:  The  range  of  that  model  term. 

We  can  draw  some  general  conclusions  from  these  quantities,  but  final  judgement 
should  be  carried  out  simultaneously  with  the  graphical  representations  and  with  the 
residual  analysis  material.  The  terms  to  be  concerned  about  are  the  large  parameters. 
If  we  included  these  in  the  overlap,  then  they  will  dominate  the  covariance  matrix 
and  lead  to  a very  unstable  solution.  If  they  are  only  varying  within  their  error  then 
they  can  be  constrained  and  as  such  will  not  effect  the  covariance  matrix  as  strongly. 
If  they  vary  too  wildly  then  it  may  be  that,  with  our  data,  particular  term  cannot  be 
realistically  modelled. 

The  terms  9,  10,  11,  12  and  9,  10  respectively  in  the  1956  and  1992  material 
vary  widely.  These  model  the  magnitude  distortion  terms:  (m  - (m)X(^2  ^7/2^ 

and  (m  - (m))r/(^2  _|_  These  terms  being  cubic  in  the  standard  coordinates  are 
very  small,  and  their  parameters  have  to  be  corresponding  large,  but  this  does  cause 
problems  for  the  inversion  process.  Since  this  term  is  a very  system  dependent  term, 
e.g.  dependent  on  the  optics  more  than  the  observing  / guiding,  it  may  be  physically 
realistic  to  constrain  it  to  one  parameter  for  all  the  plates. 

Term  13  in  the  1992  material  also  seems  to  be  very  discordant,  it  has  a small 
value,  a reasonably  small  error  and  error  standard  deviation;  yet  it  has  a large  value 
standard  deviation  and  very  large  range.  This  may  have  some  discordant  values  that  are 
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skewing  the  overall  distribution.  The  graphical  examination  also  provides  information 
with  which  to  answer  some  of  these  concerns. 

Graphical  Analysis  of  Parameters 

Figures  23,  24,  and  25  show  the  variation  of  the  first  eight  extra  parameters  for  a 
the  same  sample  of  plates  from  each  epoch.  This  is  a graphical  representation  of  the 
parameters  just  discussed.  The  y axis  represents  the  value  of  the  parameter,  and  the 
X axis  is  a running  plate  number.  The  order  of  the  plates  corresponds  to  the  date  of 
exposure.  On  each  parameter  point,  the  length  of  the  vertical  line  represents  the  formal 
error,  e.g.  the  square  root  of  the  corresponding  diagonal  term  on  the  covariance  matrix. 
The  dotted  line  represents  a two  standard  deviation  level  for  each  parameter.  The  line 
through  the  middle  is  a best  fit  horizontal  line. 

As  one  can  see  the  size,  range  and,  scatter  varies  quite  a lot.  This  ‘scatter’  is 
quantified  by  the  ratio  of  the  parameter  standard  deviation  to  its  mean  covariance  in 
table  5.  This  ratio  will  indicate  if  a parameter  is  varying  outside  its  formal  error.  If  this 
is  not  the  case  then,  providing  the  variation  is  due  to  random  causes,  we  can  assume  the 
variation  is  within  the  noise  and  should  either  be  constrained  or  dropped  completely. 
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Table  5:  Parameter  Statistics 


# 

Mean 

a 

Range 

Mean  t 

<7e 

Epoch  1900 

7 

-0 . 00155 

0 . 00219 

0 . 0104 

0 . 000787 

0.000285 

8 

-0 . 00210 

0.00173 

0 . 00869 

0.000837 

0.000318 

9 

0.148 

15.0 

58.0 

7 . 03 

2.27 

10 

11.3 

13 . 7 

58 . 0 

7.21 

2.33 

11 

0.000276 

0.000755 

0.00368 

0 . 000447 

0.000335 

12 

0.0285 

0.141 

0.788 

0 .0615 

0.0499 

13 

-28.2 

1 .21e+03 

6. 07e+03 

576. 

137  . 

14 

-490. 

1 .27e+03 

4 . 97e+03 

572  . 

135. 

15 

4.58 

10 . 7 

42.1 

5.11 

1 . 18 

16 

CM 

00 

1 

11 . 4 

53.6 

Epoch  1956 

5.15 

1 . 18 

7 

-7 . 17e-05 

0 . 00512 

0 . 0522 

0.000778 

0.000432 

8 

0 . 000473 

0.00427 

0.0395 

0 . 000783 

0 . 000386 

9 

1 . 07e+04 

3.21e+04 

3.47e+05 

6. 62e+03 

6.286+03 

10 

-562. 

3. 40e+04 

3.59e+05 

6 . 286+03 

6 . 456+03 

11 

3.72e+03 

5 . 49e+04 

6. 18e+05 

1 . lOe+04 

1.596+04 

12 

2.24e+04 

6 . lle+04 

6.51e+05 

1 .146+04 

1 . 626+04 

13 

0 . 0307 

0.710 

6.47 

0 . 146 

0 . 102 

14 

-0 . 0540 

1 . 11 

9.84 

0.220 

0.205 

15 

-1.56 

80.1 

651 . 

15.5 

15.9 

Epoch  1992 

7 

-0 . 000745 

0.0106 

0 . 112 

0 . 000666 

0.000270 

8 

0.000319 

0.00816 

0 . 0797 

0 . 000665 

0 . 000271 

9 

1 . 44e+04 

1 .21e+05 

2 . Ole+06 

4.256+03 

4.506+03 

10 

1 . 07e+04 

1 .22e+05 

2 . 12e+06 

3.986+03 

4 . 156+03 

11 

-0 . 183 

2 .74 

27 . 2 

0 . 146 

0 .113 

12 

-0.202 

3.35 

33.5 

0.213 

0.195 

13 

21.2 

302. 

4 . 21e+03 

13.2 

12.1 
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Table  5:  (Continued)  Parameter  Statistics 


# 

Mean  e / <r 

Mean  / Mean  e erg  / cr 

Term  Mean 

Term  Range 

Epoch  1900 

7 

0.358 

1 . 97 

0 . 130 

-0 . 0186 

0 . 125 

8 

0 . 485 

2.51 

0 .184 

-0 . 0252 

0 . 104 

9 

0 .468 

0 . 0210 

0 . 151 

0 . 000214 

0 . 0838 

10 

0.526 

1.57 

0 . 170 

0 . 0163 

0.0837 

11 

0.592 

0.617 

0 .444 

0 . 0397 

0.529 

12 

0.437 

0.463 

0 . 355 

0 . 0130 

0.359 

13 

0 . 478 

0 . 0489 

0 .113 

-0.00309 

0.666 

14 

0 . 450 

0.857 

0 . 106 

-0.0538 

0.545 

15 

0 . 478 

0 .897 

0 . Ill 

0 . 0132 

0 . 122 

16 

0 . 451 

1 . 60 

0 . 103 

Epoch  1956 

-0 . 0238 

0 . 155 

7 

0 . 152 

0 .0922 

0.0844 

-0 . 000861 

0 . 626 

8 

0 .184 

0.605 

0.0906 

0.00568 

0 . 474 

9 

0.206 

1 . 61 

0 .195 

1 . 14 

36.9 

10 

0 . 185 

0.0895 

0 .190 

-0 . 0598 

38.3 

11 

0.201 

0.338 

0.289 

0.286 

47.5 

12 

0.186 

1.96 

0.265 

1.72 

50 . 1 

13 

0.206 

0.210 

0 .143 

0 . 0796 

16.8 

14 

0 .199 

0.245 

0 . 185 

-0.101 

18 . 4 

15 

0 .194 

0.101 

0 .199 

Epoch  1992 

-0 . Ill 

46.2 

7 

0 . 0631 

1 . 12 

0 . 0256 

-0.00894 

1.34 

8 

0.0814 

0.480 

0 . 0332 

0 . 00382 

0.957 

9 

0 . 0351 

3.39 

0 . 0371 

1.54 

214  . 

10 

0 . 0326 

2 . 68 

0 . 0340 

1 .14 

226. 

11 

0.0532 

1.26 

0 . 0412 

-0 . 475 

70.5 

12 

0 .0634 

0 . 949 

0 . 0581 

-0.378 

62.7 

13 

0.0438 

1.60 

0 . 0401 

1.51 

299. 
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There  are  circumstances  when  that  will  not  be  the  case.  For  example  in  the  case 
of  the  straight  magnitude  parameters,  e and  f,  these  are  primarily  terms  that  result  from 
bad  guiding.  If  they  are  included  (e.g.  if  they  are  significant)  then,  because  they  are 
very  plate  / observer  dependent,  constraining  them  would  not  make  physical  sense.  The 
two  terms  may  not  however  be  significant,  if  they  vary  around  a mean  of  zero,  and 
their  standard  deviation  is  less  than  their  formal  error  they  can  be  dropped. 

The  graphs  confirm  that  the  parameters  for  magnitude  distortion  can  defiantly  be 
constrained  for  the  1992  epoch  (parameters  9 and  10),  and  probably  constrained  for  the 
1956  epoch  (parameters  9,10,11,12).  More  analysis  of  this  nature  must  be  carried  out 
to  make  a definite  decision. 

An  examination  of  this  nature  to  find  the  model  is  very  convoluted.  This  discussion 
has  shown  how  we  can  use  the  residuals  and  parameters  to  indicate  which  model  to  use. 
The  whole  choice  is,  as  we  have  said,  not  black  and  white.  After  iterating  through  tests 
on  parameters  and  residuals  simultaneously  we  will  eventually  see  a pattern  of  terms 
that  need  to  be  included.  At  some  point  the  decision  will  be  made  that  the  potential 
improvement  in  precision  does  not  warrant  further  work  on  the  model.  There  is  one 
last  reduction  improvement  we  can  make  after  this  point  to  improve  the  solution,  this 
involves  analyzing  the  output  variances. 

Below  are  the  models  we  have  decided  upon  for  this  investigation,  for  all  three 
epochs  we  included  the  basic  eight  constant  model.  For  epoch  1900  we  included  as 
extra  terms  the  following: 

{in  - {m)f  (m  - {m))r] 


(5.8) 
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Figure  23:  Plate  Parameter  Variation  for  Epoch  1900. 
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Figure  24:  Plate  Parameter  Variation  for  Epoch  1956. 

Porom  if  8 


0.015 

0.010 

0.005 

0.000 

-0.005 

-0.010 

-0.015 


T 


. Iff  •*:  '?fl 

■ ^ - +1 


■+X  -fc 


1 


+ 


50 


Porom  I 10 


112 


Figure  25:  Plate  Parameter  Variation  for  Epoch  1992. 
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with  constrained  terms  in  the  following: 


r + 0 

0 + 


(5.9) 


For  epoch  1956, 


(m  — 0 

0 (m  — {m))^T] 


(5.10) 


with  constrained  terms  in  the  following: 


(m  - (m})^(^2  4-  r]^) 
0 


0 

{m  - {Tn))r]{^^  + T}^) 


(5.11) 


Finally  for  epoch  1992, 


(m  — (m))^^  0 

0 (m  — {m))^  T] 


(5.12) 


with  constrained  terms  in  the  following: 


(m  - (m))^  (^2  ^ ^2)  _ (m))j7(^2  ^2^ 

('«- (m))7y(^^  + r?^)  (m-(m))^(^2  + ,,2) 


(5.13) 
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Overlap  Variances 

After  the  overlap  has  been  carried  out  the  stellar  positions  will  have  a error  from 
equation  3.56: 

€={sBl3  + Ea-d).  (5.14) 

The  square  of  these  quantities  will  provide  the  variance  of  the  star’s  position.  This 
variance  will  still  reflect  inherent  problems  in  finding  stellar  positions.  As  we  have 
discussed  the  positional  precision  of  a star  is  related  to  its  magnitude  and  to  its  proximity 
to  the  center  of  the  plate.  The  positions  of  the  stars  at  the  extremes  in  magnitude 
(brightest,  faintest)  and  plate  position  (in  the  corners  and  edges)  will  not  be  determined 
as  precisely  as  those  in  the  middle  of  the  range  / plate.  However,  after  one  run  we  can 
find  how  the  variance  is  correlated  with  the  stellar  parameters,  position  and  magnitude. 

In  figure  26  and  27  we  plot  the  variances  of  a nine-constant  fit,  that  includes  both  first 
order  magnitude,  and  origin  distance  for  the  AC  material.  The  variances  are  plotted  vs 
first  second  and  third  order  magnitude,  coordinate  and  origin  distance.  We  see  positive 
correlations,  and  can  fit  a least  square  polynomial  to  all  the  system  variances  and  then 
re-run  the  overlap  including  the  predicted  variance  rather  than  the  constant  variance 
we  use  in  the  first  run.  It  is  important  to  return  to  the  original  estimates  of  the  stellar 
equatorial  positions;  otherwise  the  variance  information  is  not  ‘new’  information.  If  we 
do  not  return  to  the  original  estimates  the  re-computation  will  just  reconfirm  that  our 
input  has  the  predicted  variance  variation. 

This  ends  the  discussion  on  choosing  the  model.  In  all  solutions  the  covariance 
matrix,  correlation  matrix,  matrix  efficiency,  and  condition  number  are  calculated  and 
checked  to  make  sure  the  system  is  still  stable. 
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Figure  26:  Nine-Constant  Overlap  Variances 
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Figure  27:  Nine-Constant  Overlap  Variances 
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CHAPTER  6 
PROPER  MOTIONS 


Finding  Proper  Motions 

The  proper  motion  of  a star  is  defined  as  the  time  derivative  of  a star’s  equatorial 
coordinate  in  an  (ideally)  inertial  frame.  In  this  study  the  nearly  inertial  frame  is 
provided  by  the  Astrographic  Catalogue  of  Reference  Stars  using  the  B1950  (Besselian 
Date  1950)  positions.  This  means  all  three  epochs  are,  within  the  accuracy  of  the 
reference  material,  in  the  FK4  approximation  to  an  inertial  reference  frame. 

To  keep  the  plate  parameters  and  centers  as  close  as  possible  to  those  published 
when  the  plates  were  exposed  we  precessed,  using  lAU  constants,  all  the  reference 
material  (and  by  default  the  inertial  frame)  to  1900  for  the  first  epoch.  After  we  found 
first  epoch  positions  we  then  precessed  them  back  to  B1950,  the  epoch  of  the  reference 
catalogue  and  the  second  and  third  observations. 

The  difference  in  the  stellar  positions,  after  this  procedure,  between  the  three  epochs 
is  due  to  the  proper  motion  of  the  respective  stars.  The  relationship  of  positions  (a, 6) 
at  any  given  date  (t),  with  the  stellar  positions  (oo,  <5o)  and  the  proper  motions  (/t)  at 
the  date  of  observation  (to),  are 

a = Oo  + Ha(t  -to)  6 = 6o  + fj.6(t -to)-  (6.1) 

Note  that  the  proper  motion,  ft,  has  no  time  suffix;  this  is  because  we  assume  uniform 
motion  and  therefore  — fig- 

To  obtain  the  best  estimates  for  the  position  and  proper  motion  we  should  find 
these  quantities  at  the  position  variance  weighted  epoch.  At  this  point  the  position  and 
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proper  motion  estimates  are  uncoupled  and  therefore  uncorrelated.  In  this  example  the 
three  position  estimates  have  a variance  ratio  of  approximately  4:1:1  in  1900:1956:1992 
respectively.  The  variance  weighted  average  epoch  is  then  (1900+(4*1956)+(4*1992))/9 
« 1965. 

However,  the  plates  were  not  all  exposed  at  the  same  time.  It  has  already  been 
pointed  out  that  the  AC  material  was  exposed  over  a period  of  10  years  and  the  second 
epoch  McCormick  observations  took  over  a year.  In  the  overlapping  plate  method  we 
assume  that  a star  on  different  plates  (within  the  same  epoch)  has  the  same  position. 
However,  because  of  the  time  difference  between  exposures  this  will  not  be  the  case. 

One  way  to  alleviate  the  problems  the  above  assumption  could  cause  is  the  follow- 
ing; after  the  overlap,  rather  than  using  the  /3  value  to  correct  the  input  and  get  one 
mean  weighted  position,  actually  use  the  plate  parameters  to  find  the  position  of  a star 
on  each  plate.  The  strength  of  the  overlap  enters  into  this  estimate  of  the  star’s  position 
via  the  plate  parameters.  However,  the  actual  star  position  is  no  longer  confined  to 
just  the  average  for  that  epoch,  it  can  rather  assume  the  value  derived  from  each  plate 
along  with  the  exposure  date  of  that  plate.  So  if  two  plates  that  overlap  have  a date  of 
exposure  difference  of  ten  years,  then  the  star’s  position  will  be  allowed  to  reflect  this. 

There  is  a more  rigorous  method.  Eichhorn  (1993)  has  derived  a method  that  allows 
complete  generality  in  the  position  by  introducing  the  proper  motion  as  an  additional 
parameter  at  the  beginning  of  the  reduction.  In  this  reduction  the  time  difference  is 
most  acute  in  the  AC  material,  which  is  also  the  material  with  the  lowest  positional 
precision.  The  increased  computational  effort  required  to  incorporate  Eichhorn’s  new 
method  is  therefore  unwarranted  and  the  above  approximate  approach  is  sufficient. 
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Therefore  we  find  the  proper  motion  from  a least  squares  analysis  of  all  occurrences 
of  a star  in  every  exposure.  Following  the  discussion  in  chapter  two  we  first  set  up 


the  condition  equations, 

/ a - (ai  + /ia(1965  - ti))  \ 
6-(^i  + /ii(1965-<i)) 
a - («2  + ^a(1965  - <2)) 

(5  - ((?2  + /ii(1965  - <2)) 


F = 


a - (««  + /^a(1965  -<„)) 

V 5-(5„  + /i^(1965-<„))  / 


/ a - ai  + Haiti  - 1965)  \ 
6 — 61  -{-  Hiiti  ~ 1965) 
a - 02+  Ha{t2  - 1965) 

S — 62  + Hi{t2  ~ 1965) 

a - «„  + Ha{tn  - 1965) 

V (5  - 6„  + Hs{tn  - 1965)  J 


(6.2) 


where  &nd  fis  are  the  target  quantities,  ai, 02,... a„  and  ^2, •••^2  are  the 

‘observations’  from  the  plate  solutions  at  dates  <1,^2,  •••(3.  respectively. 

Following  equations  2.12  and  2.13,  we  define  the  two  quantities:  A = 

i%)x=x  a=a  and  X = (^)j._j.  then  the  corrections  to  the  parameters  a (equation 
2.20)  are 

a = - a'^(X(tX'^)"V  ”\'^(X(tX'^)~^Fo.  (6.3) 

As  in  the  overlap  we  have  the  simplification  that  each  equation  contains  only  one 
observation  and  therefore  X = - I.  The  matrix  A is 


A = 


d¥ 


d{a,b,Hoc,H6) 


/I 

0 

T\ 

0 \ 

0 

1 

0 

n 

1 

0 

T2 

0 

0 

1 

0 

T2 

1 

0 

Tn 

0 

\0 

1 

0 

Tn) 

(6.4) 


where  t is  the  time  difference  t,-1965.  Assuming  a first  approximation  of  zero  for  the 
four  unknowns,  then  Fo  is 


F„  = 


-02 

-^2 


V 

-On 


(6.5) 
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where  t is  the  time  difference  t,-1965.  Assuming  a first  approximation  of  zero  for  the 
four  unknowns,  then  Fq  is 


F„  = 


-tt2 


—Oin 


(6.5) 


The  parameters  simplify  to  a = 

= -[aV 

-A]- 

The  first  term 

simplifies  to 

0 

[aT<t-'a]  = 

0 

E<^J 

^ 0 

(6.6) 

0 

0 

V 0 

E^r 

0 

where  all  terms  are  summed 

over  all  observations.  This 

is  a (4  X 4) 

matrix. 

The 

second  term  simplifies  to  the  (4x1)  vector: 

( \ 


V ) 

The  parameter  errors  will  be 

(diag[A'^(T“^A]  ’ 

and  the  observation  errors  (residuals)  will  be  £ = (Fq  — Aa). 


(6.7) 


(6.8) 


Results 

Table  (9)  lists  the  calculated  position  and  proper  motion  for  the  epoch  1965  equinox 
1950B.  All  stars  that  were  found  in  at  least  two  epochs  are  included.  Of  the  original 
2600  stars  from  the  Astrographic  Catalogue  Plates  1678  were  found  in  at  least  one 
more  epoch.  The  others  are  predominantly  in  the  edge  regions  not  covered  in  the 
McCormick  Observations. 


121 


From  the  first  epoch  there  were  6285  observations,  from  the  second  epoch  there  were 
35,378  and  from  the  last  epoch  33,352.  The  positional  variances  were  approximately 
in  the  ratio  4:1:1  for  the  1900:1956:1992  epochs  respectively.  This  is  expected  if  one 
considers  the  higher  measurement  precision  and  longer  focal  length  in  the  second  and 
third  epochs. 

Table  (6)  lists  the  mean  error  and  standard  deviation  of  the  formal  errors  in  the 
final  solution. 


Table  6:  Final  Solution  Statistics 


Quantity 

Right  Ascension 

Declination 

Position  Error  Mean,  Standard 
Deviation.  " 

0.0142,  0.013 

0.050,  0.046 

Proper  Motion  Error  Mean, 
Standard  Deviation,  "/century 

0.0058,  0.0035 

0.0217,  0.0139 

One  way  to  internally  check  the  result  is  to  see  how  the  results  of  comparing 
the  1900  material  to  only  one  of  the  other  two  epochs  differs  from  that  of  the  least 
squares  answer.  Figure  28  presents  the  results  of  a very  simple  calculation  that  takes 
the  difference  in  the  stellar  position  at  the  two  epochs  and  extends  it  out  by  a factor 
of  200.  The  star  is  at  the  tail  of  the  proper  motion.  As  one  can  see  if  one  takes  into 
account  the  fact  that  the  epoch  difference  is  almost  twice  as  much  for  the  1900-1992 
than  the  1900-1956  observations,  these  graphs  are  qualitatively  the  same. 


declination/d  declination/d 
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Figure  28:  Proper  Motions  Using  Just  Two  Epochs 
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Figure  29  shows  the  least  squares  solution  for  proper  motions  over  a 10,000  year 
period.  These  also  reflect  the  general  pattern  seen  in  the  last  two  graphs.  The  high 
proper  motion  star  in  the  middle  of  the  picture  is  HD  36443  a 8.6  magnitude  G5  star 
that  is  not  believed  to  be  part  of  the  association. 

Note  that  none  of  the  belt  stars  are  included  because  they  were  not  measured 
in  the  Astrographic  Catalogue.  The  positions  for  stars  as  bright  as  these  are  very 
difficult  to  determine  from  photographic  plates  because  of  the  sizes  of  their  images. 
For  informational  purposes,  their  positions  are  marked  by  their  names:  e,  S , and  C- 
In  the  second  graph  in  figure  29  we  show  the  proper  motion  diagram  of  our  sample. 
If  the  association  members  had  a significant  common  proper  motion,  then  they  would 
group  together  in  the  same  portion  of  the  diagram.  The  mean  proper  motion  of  this 
group  would  then  define  the  mean  proper  motion  of  the  association.  However,  an 
examination  of  the  diagram  shows  no  such  structure,  and  the  membership  is  not  easily 
obtainable  from  just  the  proper  motions.  After  further  analysis  of  the  proper  motion 
components,  we  felt  that  membership  determination  using  proper  motions  would  require 
an  examination  beyond  the  scope  of  this  dissertation. 
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Figure  29:  Proper  Motion  Plots  Final  Result 


125 


A good  way  to  judge  the  accuracy  of  our  results  is  to  compare  the  positions  and 
proper  motions  with  those  of  the  Astrographic  Catalogue  of  Reference  Stars.  This  is  still 
an  internal  check  because  these  reference  stars  were  used  as  our  original  input.  It  may 
however  highlight  systematics.  Figure  30  shows  the  proper  motions  of  the  reference 
stars  in  our  region  taken  directly  from  the  ACRS.  Again  there  is  good  qualitative 
agreement  between  this  figure  and  the  first  graph  in  figure  29. 


Reference  Star  Proper  Motions 


Figure  30;  Proper  Motions  of  Reference  Stars  in  our  Sample 
By  finding  the  difference  between  the  reference  stars  positions  and  proper  motions 
in  1965  and  our  calculated  positions,  we  can  find  the  ‘residuals’  of  our  results.  This 
can  also  be  done  by  comparing  our  results  with  AGK3U,  a updated  version  of  the 
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AGK3  catalogue.  We  have  341  stars  overlapping  with  the  ACRS  and  206  stars  with 
the  AGK3U.  Figure  31  shows  the  histograms  of  the  difference  between  our  positions  / 
proper  motions  and  these  two  catalogues.  The  curve  represents  a best  fit  gaussian. 

The  striking  feature  of  these  plots  is  that  the  residual  in  declination  is  distributed 
over  two  times  the  range  of  the  right  ascension  residual.  There  could  be  a number  of 
intrinsic  reasons  for  this  which  may  be  related  to  the  actual  observations  and  measuring. 
For  example  the  observations  were  all  hand  guided,  and  the  guiding  corrections  made 
during  an  actual  exposure  were  nearly  always  in  right  ascension.  The  McCormick 
telescope  is  equatorial  mounted  telescope,  thus  if  the  telescope  had  an  axes  alignment 
problem,  then  there  would  be  a drift  in  declination.  This  may  be  imperceptible  to  the 
observer,  but  would  lead  to  a systematic  guiding  error  in  the  declination  direction  and 
hence  a larger  variance  in  that  coordinate.  Another  possibility  for  the  larger  declination 
residual  may  be  due  to  the  measuring  machine.  The  plates  were  all  placed  on  the 
machine  so  that  the  x and  y coordinates  always  had  the  same  orientation  with  respect 
to  the  measuring  axes.  If  there  was  a continual  slip  in  the  y direction  this  would  also 
cause  a large  variance  in  the  declination  because  this  axis  during  observation  would 
be  aligned  with  that  coordinate.  A re-analysis  of  both  observing  device  / procedure 
and  the  measuring  device  / procedure  is  needed  to  distinguish  between  these  and  other 
possibilities. 
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Figure  31:  Differences  between  these  Results  and  External  Catalogues 
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This  ends  our  discussion  of  determining  the  proper  motions  in  Orion.  The  re- 
sults, data  and  programs  will  be  available  by  anonymous  file-transfer-process  from 
astro.ufl.edu.  Next  we  will  discuss  the  use  of  these  proper  motions  and  the  options 
for  future  work. 

Using  the  Proper  Motions 


Historically  proper  motions  of  associations  and  clusters  have  been  used  to  find  the 
expansion  age  and  kinematical  properties  of  the  association  / cluster.  These  proper 
motions  need  to  be  included  with  spectroscopic  and  radial  velocity  data  to  establish 
membership  criteria.  Once  the  members  are  determined,  the  expansion  age  and  other 
astrophysically  interesting  properties  such  as  the  H-R  diagram,  initial  mass  function, 
and  kinematical  properties  can  be  found. 

As  stated  before,  it  is  beyond  the  scope  of  this  investigation  to  examine  the 
membership  criteria,  but  we  shall  use  our  data  to  find  what  the  simple  linear  expansion 
age  would  be  assuming  the  association  is  expanding.  We  have  adopted  the  method 
developed  by  Blaauw  (1952)  to  determine  a linear  expansion.  He  showed  that  should 
an  expansion  exist  and  provided  that;  1)  the  present  dimensions  of  the  group  are  larger 
than  the  original  dimensions  and  2)  all  the  motions  in  the  group  started  at  the  same  time 
and  are  uniform,  then  there  should  be  a linear  relationship  between  the  component  of  the 
proper  motion  and  the  corresponding  coordinate.  One  must  also  allow  for  perspective 
contraction  or  expansion  caused  by  the  relative  motion  of  the  Sun. 

Using  this  method  Lesh  (1968)  found  an  expansion  age  of  4.5  billion  years  for 
the  northwest  region  of  I-Orion.  This  examination  was  based  on  16  stars  which  were 
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not  found  in  our  sample.  However  there  have  been  two  other  studies  (Warren  and 
Hessler,  1977  and  Gieseking,  1983)  that  have  found  possible  members  which  are  within 
our  data  set.  We  will  adopt  their  membership  choices  and  use  our  proper  motion 
components  to  investigate  the  expansion.  The  Warren  and  Hessler  study  used  primarily 
a photometric  determination  of  membership,  while  the  Gieseking  study  was  a radial 
velocity  investigation  of  highly  probably  members.  We  will  consider  these  as  studies 
of  the  Orion  members  using  spectroscopic  and  radial  velocity  criteria,  respectively. 

Spectroscopic  Criteria 

There  is  a sizable  amount  of  spectroscopic  information  available,  however,  most 
of  it  is  concerned  with  the  bright  stars.  The  above  mentioned  study  by  Warren  and 
Hessler  in  1977  (hereafter  WH)  collected  all  the  spectroscopic  criteria  available  and 
combined  this  with  ubvy  observations  to  produce  a membership  ranking  for  504  stars 
in  the  association,  of  which  107  overlap  with  this  study.  The  stars  are  all  ranked 
on  their  probability  of  membership  in  radial  velocity,  proper  motion  and  photometric 
properties  as  a,  b or  c.  An  ‘a’  represented  a high  probability  of  membership,  while  a 
‘c’  represented  a low  probability. 

In  figure  30  we  plot  the  proper  motion  components  vs.  their  respective  coordinates. 
A variance  weighted  least  squares  fit  to  the  data  is  shown.  The  first  two  graphs  represent 
very  high  probability  stars  which  were  classified  as  ‘a’  by  WHs’  membership  criteria  in 
all  three  respects;  proper  motion,  radial  velocity,  and  photometry.  If  the  straight  line  fit 
has  a positive  slope,  then  this  indicates  an  expansion;  while  a negative  slope  indicates  a 
contraction.  As  one  can  see,  the  slopes  are  not  strongly  positive,  they  are  either  basically 
flat  or  negative.  The  slopes  and  the  standard  deviations  of  the  least  squares  fit  are  shown 
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in  table  7.  None  of  the  slopes  are  above  the  standard  deviation,  so  their  significance  must 
be  questioned.  However,  should  these  slopes  actually  signify  true  kinematic  properties 
of  the  association,  then  in  actual  fact  the  association  is  contracting. 

It  must  be  noted  that  this  is  a very  simplified  analysis  that  can  be  questioned  on 
many  different  grounds.  For  example  the  membership  criterion  was  adopted  wholesale 
from  WH;  and  by  using  more  recent  observations,  better  membership  criteria  may  be 
possible.  We  have  also  not  allowed  for  a distance  gradient  as  proposed  by  WH,  because 
there  was  not  enough  available  distance  information  for  this  small  sample  size  to  make 
this  consideration  useful.  Also,  there  is  a possibility  that  galactic  rotation  might  distort 
the  proper  motions.  All  of  these  considerations  were  considered  outside  the  scope  of 
this  investigation  but  will  be  interesting  to  consider  as  future  projects. 


Table  7:  Lxast  Squares  Slopes  using  WH  Criteria. 


Coordinate 

Number  of 
stars  in  sample 

Slope  deg'^ 
year'^ 

Standard 

Deviation 

Right  Ascension,  all  criteria  =a. 

26 

-0.0049 

0.8881 

Declination,  all  criteria  =a. 

26 

-0.2581 

0.7999 

Right  Ascension,  photometric 
criteria  =a. 

103 

-0.0032 

1.1089 

Declination,  photometric 
criteria  =a. 

103 

-0.0598 

0.7977 
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u vs.  a Best  Volues 


Figure  32:  Linear  Relations  between  yua  and  a and  and  S using  WH  Membership  Criteria 
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Radial  Velocity  Criteria 


Radial  velocities  can  also  be  used  to  determine  association  membership.  However, 
measuring  radial  velocities  is  very  time  consuming  and  to  date  only  a handful  of  the 
stars  in  this  study  had  their  radial  velocities  measured.  Due  to  the  inherent  observational 
difficulty  of  measuring  radial  velocities,  these  are  also  predominantly  bright  stars.  If 
radial  velocities  from  different  observatories  are  used,  then  an  investigation  of  individual 
observatory  corrections  would  be  required.  Therefore,  we  will  adopt  the  results  of 
Gieseking’s  (1983)  study  of  66  of  high  probability  Orion  Ib  members.  Gieseking  felt 
there  was  a bimodal  distribution  in  the  group’s  radial  velocity  distribution  centered  on 
4.5km/s  (relative  velocity). 

Of  the  66  stars  in  his  study  50  are  also  in  this  study.  The  expansion  was 
again  investigated  by  plotting  the  proper  motion  components  against  their  respective 
coordinates,  as  shown  in  figure  33.  In  these  graphs  the  symbols  *•’  represent  the  stars 
with  a velocity  below  4.5  km/s  and  ‘+’  those  greater  than  4.5  km/s.  The  three  lines 
represent  a least  squares  best  line  fit  to  the  stars  with  velocities  less  than  4.5  km/s 
(dotted  line),  those  with  velocities  greater  than  4.5  km/s  (dashed  line),  and  all  the  stars 
in  the  sample  (solid  line).  The  dashed  circles  represent  the  variance  of  each  star’s  proper 
motion  measurement.  The  slopes  and  the  standard  deviations  of  the  least  squares  fit 
are  shown  in  table  8. 

Again  nothing  conclusive  can  be  stated  because  the  standard  deviation  are  greater 
than  the  slope  of  the  positions.  However  there  is  a consistent  expansion  in  both 
coordinates  and  it  is  possible  to  calculate  an  expansion  using  these  slopes. 

Because  Orion  is  receding  almost  directly  in  the  direction  of  the  solar  antapex  we 
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can  assume  that  the  recession  velocity  will  give  an  apparent  contraction  (Blaauw  1952, 
Lesh  1968)  of 


7T  V,  1 
180  4.74  r 


(6.9) 


where  Vr  is  the  velocity  of  recession  and  r is  the  distance.  Using  a recession  velocity 
of  20.8  kms~^  and  a distance  of  400  pc  we  arrive  at  a contraction  due  to  recession  of 
0.019  degrees  / century.  Using  a weight  least  square  slope  of  all  the  observations  in 
both  coordinates  and  adding  in  this  recessional  term  we  find  an  expansion  age  of  1.22 
X 10^  years  with  a standard  deviation  of  200%.  This  is  basically  consistent  with  what 
one  would  expect  but  in  fact  totally  inconclusive. 


Table  8:  Least  Squares  Slopes  using  Gieseking  Criteria. 


Coordinate 

Number  of  stars  in 
sample 

Slope  deg'^  year'^ 

Standard  Deviation 

Right  Ascension, 
all  criteria  =a. 

55 

0.468 

0.936 

Declination, 
all  criteria  =a. 

55 

0.145 

0.756 

Average  of  both 
slopes 

55 

0.277 

0.582 

'/century  /x  "/century 
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Figure  33:  Linear  Relations  between  ficx  and  a and 
and  6 using  Gieseking  Membership  Criteria 
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It  appears  that  using  the  best  currently  published  membership  criteria  and  a simpli- 
fied analysis  that  no  conclusive  proof  of  either  expansion  or  contraction  can  be  obtained. 
This  will  be  discussed  again  in  the  section  on  future  work. 

Future  Work 

The  proper  motions  estimates  derived  in  this  investigation  may  be  improved  with 
more  information,  primarily  color  and  additional  observations.  The  models  we  checked 
never  included  color  terms  because  we  did  not  have  this  information  for  the  field  stars. 
When  the  Guide  Star  Catalogue  II  is  available  this  will  change  as  this  catalogue  is 
complete  down  to  12th  magnitude,  in  combination  with  the  first  catalogues  will  give 
two  color  magnitudes.  This  is  obviously  not  as  useful  as  spectral  types  but  will  suffice 
for  testing  the  material  for  any  color  correlations.  In  addition  to  providing  colors,  the 
two  guide  star  catalogues  will  provide  two  extra  epoch  observations.  The  estimates  will 
not  be  as  accurate  because  the  guide  star  has  a high  error  budget  and  used  astrometrically 
inferior  Schmidt  plates,  but  they  will  still  provide  valid  estimates  of  the  stellar  positions. 

The  reduction  procedure  used  in  this  investigation  can  be  improved  upon.  Some 
of  the  parameters,  notably  the  focal  length,  are  very  strong  candidates  for  stochastic 
constraints.  This  will  allow  for  less  variance  in  the  final  solution.  We  could  also  include 
proper  motions  from  step  one.  Rather  than  reducing  the  three  epochs  separately  one 
can  reduce  them  simultaneously,  the  proper  motions  will  then  be  parameters  from  that 
reduction.  This  will  provide  for  more  internally  consistent  estimates  of  all  the  stellar 
parameters.  Finally,  the  model  we  have  chosen  can  probably  be  improved  with  more 
investigation.  This  is  always  an  option  but  the  expected  gains  are  less  than  those  from 
incorporating  one  of  the  improvements  above. 
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This  study  provides  consistent  proper  motions  estimated  for  a large  sample  of  stars 
in  the  region  of  the  Orion  association.  These  need  to  be  investigated  in  conjunction 
with  all  the  currently  available  spectroscopic  and  radial  velocity  criteria  to  determine  as 
large  a membership  sample  as  possible  which  is  free  from  contamination.  This  sample 
should  then  be  examined  with  the  methods  outlined  above,  using  space  velocities,  and 
using  proper  motions  in  galactic  coordinates.  These  examinations  should  be  applied  to 
the  association  as  a whole  and  to  the  respective  subgroups  individually.  This  study  has 
shown  that  when  considering  previously  determined  member,  nothing  conclusive  can 
be  drawn  about  the  expansion  or  contraction  of  the  association. 
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Table  of  Results 

Table  9 lists  the  proper  motions  and  positions  for  stars  found  in  the  Astrographic 
Catalogue  and  at  least  one  of  the  McCormick  Epochs.  The  epoch  of  the  position  is 
1965  and  the  equinox  is  1950B.  Listed  are: 

1.  AC  # — the  running  number  for  the  star  in  the  1900  AC  plates 

2.  ACRS  # — the  ACRS  number  if  the  star  is  also  part  of  the  reference  catalogue 

3.  N — the  number  of  images  found  for  the  star 

4.  m — the  magnitude  from  the  AC. 

5.  a — the  right  ascension  in  hours,  minutes  and  seconds 

6.  €a  — the  right  ascension  position  error  in  seconds 

7.  6 — the  declination  in  hours,  arcminutes  and  arcseconds 

8.  — the  declination  position  error  in  seconds 

9.  /xq;,  — the  right  ascension  proper  motion  in  seconds  / century 

10.  e//Q  — the  right  ascension  proper  motion  error  in  seconds  / century 

11.  — the  declination  proper  motion  in  arcseconds  / century 

12.  — the  declination  proper  motion  in  arcseconds  / century 


Table  9:  Positions  and  Proper  Motions  for  Stars  Found  in  at  least  two  Epochs. 
AC#  ACRS#  N m a ^ es  na 
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